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We prove that many cosmological models characterized by vectors nonminimally coupled to the 
curvature (such as the Turner- Widrow mechanism for the production of magnetic fields during in- 
flation, and models of vector inflation or vector curvaton) contain ghosts. The ghosts are associated 
with the longitudinal vector polarization present in these models, and are found from studying the 
^\ • sign of the eigenvalues of the kinetic matrix for the physical perturbations. Ghosts introduce two 

' main problems: (1) they make the theories ill-defined at the quantum level in the high energy/sub 

, horizon regime (and create serious problems for finding a well behaved UV completion) ; (2) they ere- 

ate an instability already at the linearized level. This happens because the eigenvalue corresponding 
to the ghost crosses zero during the cosmological evolution. At this point the linearized equations for 
^ ^ the perturbations become singular (we show that this happens for all the models mentioned above). 

. We explicitly solve the equations in the simplest cases of a vector without vev in a FRW geometry, 

and of a vector with vev plus a cosmological constant, and we show that indeed the solutions of the 
linearized equations diverge when these equations become singular. 
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I. INTRODUCTION 



. Although the WMAP measurements of the Cosmic Microwave Background (CMB) strongly support the inflationary 
' paradigm several studies pointed out some peculiar features in the data that seem at odd with the simplest 
inflationary predictions. These so called 'anomalies' include the low power in the quadrupole moment the 
O ' alignment of the lowest multipoles a ~ 5° cold spot with suppressed power |6|, an asymmetry in power between 
^ \ the northern and southern ecliptic hemispheres 0] , and broken rotational invariance Q . ^ The significance of some 
^ ■ of these effects has increased in the latest studies [1, based on the WMAP-5 years data. These observations 
I I ' have motivated a number of studies both on data analysis (how to construct estimators that can assess the degree of 
violation of statistical isotropy, see eg. [ll[) and on theory (how to construct models that reproduce these features). 
^""^ One can for instance attempt to ascribe the departure from statistical isotropy to initial conditions at the onset of 
^ inflation either on the background evolution [H, [13, [III , or on a super-horizon isocurvature mode [l^ . ^ 

Ref. [l^ associates the super-horizon mode with a long cosmic string. In general, however, in these models the 
ly-j , breaking of statistical isotropy is built-in as an initial or "boundary" condition, rather than being predicted from 
■ first principles (e.g., from a given lagrangian). A second major problem that the above models suffer is that inflation 
rapidly removes any initial background anisotropy [T^ , and blows to unobservably large scales any perturbation that 
was present at its onset. Therefore, these proposals require a minimal (and, therefore, tuned) amount of inflation. To 
' improve over these two problems, one may consider introducing in the inflationary model some nonminimal ingredients 
. that contrast the rapid isotropization caused by the inflaton. This has been realized through the addition of quadratic 
curvature invariants to the gravity action |20| , with the use of the Kalb-Ramond axion [2l| , or of vector fields [l^] ■ ^ 
The present work continues a series of previous papers in which we studied the stability of some of the models with 
^ ' vector fields. In Ref. [1^, we presented a general discussion valid for three groups of models, characterized by (i) a 
?H " potential V {A^) for the vector [l^l, (n) a fixed spatial norm of the vector, enforced by a lagrange multiplier [2^, or 
^ (iii) a nonminimal coupling of the vector to the scalar curvature [13, [H, [H, [131 ■ We showed that all these models 
have ghost instabilities. Although few explicit computations were given in the Letter [2H |. we explained the physical 
reasons behind the instability: the terms that characterize the above models break the U(l) symmetry that would 
be otherwise associated to the vector field. This introduces an additional polarization (the longitudinal vector mode) 



^ Ref. 9] found an upper limit on the anisotropy, which is compatible with the result of Q. 

^ The super-horizon mode breaks the translational, and therefore can give rise to the asymmetry claimed in Q; the Bianchi-I backgrounds 
considered in Jl2l. Il4l llSH have instead planar symmetry, and have the correct structure to explain the violation of rotational invariance 
observed by Q. General expressions for the correlation (a^m °-iim') between different multipoles were given in [12| | for generic breaking 
of rotational invariance, and in llTll for generic breaking of translational invariance. 

^ While the models studied here have been proposed in the context of primordial inflation, vector fields with nonvanishing spatial vacuum 
expectation value (vev) have been also employed as sources of the late time time acceleration |23ll . See also |24| for a different model of 
anisotropic dark energy. 
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that, in these models, turns out to be a ghost. In Ref. [31| we provided the explicit computations for the case (ii). 

In the present work, we present explicit computations for the models in the group (iii). The first of such models 
p7l | is a well known mechanism for the generation of magnetic fields. Among other things, it was pointed out there 
that, for the specific action 
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(1) 



the equations for the vector field (more precisely, those governing the vev, and the transverse polarizations) are 
the same as those of a scalar field minimally coupled to the curvature. This fact was also exploited in a series 
of recent papers, motivated by the WMAP anomalies. The work [2^ studies inflation driven by N vector fields. 
The simplest realization of 28] is characterized by three mutually orthogonal vectors with equal vev: this provides 
a homogeneous and isotropic (Friedmann-Robertson- Walker, FRW) evolution; however, one can envisage the more 
complicated situation in which a large number N of fields is present, with random orientation. This provides a nearly 

isotropic expansion, with a naturally small, O ^l/\/iV^, anisotropy. Ref. [1^ provides a simpler version of the same 

idea, in which a single vector field plus an inflaton field (in a sense, replacing the average effect of the vectors) are 
present. Finally, Refs. -iO used the vector as a curvaton field, in order to produce a nearly scale invariant spectrum 
of primordial perturbations. 

The presence of ghosts has not emerged in previous stability studies of such models, due their partial nature. 
The original work [23| studied the evolution of modes of the magnetic field, associated to the transverse photon 
polarizations, but did not discuss the role of the additional longitudinal mode. The results of 30] are based on the 
6N formalism, which computes the classical evolution of super-horizon modes, assuming that the quantum theory is 
under control (as it happens in the case of scalar field inflation, for which it was developed). Ref. [s^l studied the 
gravity waves (GW) in the model of vector inflation ^2§\, assuming that the coupling of these modes with the other 
perturbations - which, for this model, is present already at the linearized level, see eq. (I77p below - can be disregarded. 
Finally, Ref. [S^ studied the linearized equations of motion for vector inflation either in the short wavelength, or in the 
long wavelength regime. ^ The study of the linearized equations alone does not allow to see whether a perturbation 
is a positive or negative energy mode, and for this reason the presence of ghosts does not appear in this analysis. ^ 
Moreover, as we show below, the linearized equations of motion for the perturbations become singular close to horizon 
crossing, in a regime where neither the long nor the short wavelength analyses apply. In Section IVDl we provide a 
more detailed discussion of some claims made in [33] . 

It is important to stress that the ghost instability takes place also for those of the above models which have a FRW 
background. The presence of a ghost in these models is indeed due to the specific sign of the effective mass term 
A/2 — —R/6 for the vector induced by the coupling to the curvature. Indeed, as we showed in [2^, a negative mass 
squared for a vector results in a ghost, and not simply in a tachyon as in the scalar case. We stress that a massless 
vector has only the two transverse polarizations; it is therefore not surprising that the mass term controls the nature 
of the longitudinal mode. Another example on how the mass term controls the nature of a mode is given by the more 
complicated case of the graviton mass m\h^h'^ — m^hfj^^h^'^ . It is well known that, unless mi = ra2 (namely, the 
Fierz-Pauli mass term [34]) the linearized spectrum of perturbations has a ghost. 

The most obvious problem associated to a ghost is the instability of the vacuum. If a ghost is coupled to a 
normal field (and, in all the above theories, there are at least gravitational couplings), the vacuum will decay in 
ghost-nonghost excitations, with a rate which is UV divergent (since the final state quanta can have arbitrarily large 
momentum without violating energy conservation). To avoid the associated instantaneous vacuum decay, theories 
wi th g hosts are thought to be consistent only as effective theories, valid below some energy scale A. It has been shown 
in [33 that, a theory which has a ghost today coupled gravitationally to positive energy fields is phenomenologically 
viable only if A < MeV , otherwise we would see signatures of the vacuum decay in the diffuse 7— ray background. A 
stronger coupling would result in a tighter bound on A. All the above models discuss physics at much greater energy 



Moreover, we show below that also the classical evolution for the longitudinal mode leads to a divergence at super-horizon, when the 
total mass of the vector vanishes. 

® In inflationary models, the wavelength of a perturbation grows nearly exponentially, while the Hubble rate H is nearly constant. 
Therefore, at sufficiently early times, the wavelength of any mode is smaller than the horizon scale , while the opposite is true at 

sufficiently late times. The moment at which the two scales are equal is dubbed "horizon crossing"; we remark that perturbations of 
different size cross the horizon at different times (larger modes, exit the horizon earlier during inflation). 

^ Cf. the case of a scalar field, with lagrangian £ = ± (^<j)^ — m^<j)^^ /2 (dot denotes derivative with respect to time, and we ignore spatial 

dimensions); the equation of motion is ^ + <j> = for either overall sign of the lagrangian, so that no instability appears from the 
equation. However, the minus sign in the kinetic term corresponds to a negative energy field (a ghost). 
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scales than MeV, so, a necessary (but, as we will see, not sufficient) condition for them to be consistent is that the 
ghost vanishes at late times. The easiest way to achieve this is to add a positive contribution > to the mass 
term, Af^ = —R/6 + m? , and to require that m? is lower than —R/Q during inflation, but greater than it today. In 
fact, this is already built in in the models of vector inflation (28l. [29|. for which m? controls the slow roll evolution 
of the vev of the vector. The same can be also done in the case of vector curvaton, where it is found that the mass 
provides an O (rn^/H^) departure from scale invariance {H being the Hubble rate during inflation) ^o[. For the case 
of photons, there are strong upper limits on the allowed value for m (see the limits on the PDG Particle Listing ^36^). 

These limits are however less stringent than the present value of i?/6 = O (lO^'^'^eV) . 

Even if it is straightforward to eliminate the ghost (s) from the present spectrum, this does not eliminate the 
instability of these models during or after the inflationary stage. We will comment on problems that arise at the 
nonlinear level in the Discussion Section IVII We stress here that the ghost instability of these models manifests itself 
already at the linearized level. The spectrum of the theory is obtained by computing the quadratic action for the 
dynamical perturbations Yi around the given background. We do so by following, and generalizing to the present 
case, the exact same steps that are done in the standard computation of scalar fleld inflation to obtain the canonical 
modes of the system [37| . The quadratic action is characterized by a kinetic Kij Y* Yj term, a mass term, and a 
mixed term (see eq. ([Ml))- The eigenvalues of K control the nature of the physical modes (i.e., a mode is a positive 
/ negative energy excitation if the corresponding eigenvalue is positive / negative). The kinetic matrix depends on 
background quantities, and thus on time. As a consequence, its eigenvalues depend on time, and the nature of a 
mode can change during the background evolution. From our computations, we find two different types of behavior, 
according to whether the vev of the vector is or is not vanishing: 

• (A) — 0: one eigenvalue of K is negative during most of the sub-horizon regime; it changes sign at some moment 
close to horizon crossing, without passing through zero (it does so by diverging) . It crosses zero later on, when 
the total mass vanishes. 

• (A) 7^ 0: one eigenvalue of K starts positive, but crosses zero at some point close to horizon crossing, and 
remains negative for some amount of time. 

The different behaviors are due to the fact that, for a nonvanishing vev, the perturbations of the vector are mixed 
with those of the metric at the linearized level (i.e. in the quadratic action), and this affects the spectrum of the 
theory. However, we see that the mixing with gravity does not remove the ghost. The case of vanishing vev applies to 
the computations of [131 and [l^l . ^ The case of a vector with a vev applies to the model of |29|| . The vector inflation 
model of [1^ is more complicated, since it contains an arbitrary number iV of vectors (and, hence, of ghosts). We 
study the simplest realization, with three mutually orthogonal vectors. We find that, in this case, the model has three 
ghosts. Two eigenvalues behave as in the (A) = case just mentioned, while the third eigenvalue behaves as in the 
(A) case. Hence, it appears that in this case only one linear combination of the ghosts is affected by the coupling 
to gravity. 

As we show below, the system of linearized equations for the perturbations becomes singular when one of the 
eigenvalues of the kinetic matrix crosses zero. Correspondingly, we expect that the linearized solutions diverge at this 
moment. We explicitly solved the equations in the case of a single vector with no vev, and in the case of one vector 
with nonvanishing vev plus a cosmological constant. We did not solve the equations for the case of vector inflation (due 
to technical difficulties: the system contains 18 gauge invariant perturbations in its simplest realization). However 
we stress that, also for this model, we explicitly proved that (i) there are ghosts (which is by itself enough to pose 
serious doubts on any prediction obtained from this model; see the Discussion Section), and (ii) the equations for the 
perturbations become singular at some finite moments of time. 

The plan of the paper is the following. In Section |TT1 we review the basic mechanism for the models of generation 
of primordial magnetic fields, vector inflation, and vector curvaton. In Section [III[ we outline the computation of the 
quadratic action for the physical modes of the system. We also show there why vanishing eigenvalues of the kinetic 
matrix K result in singular linearized equations for the perturbations. In Section IIVI we present the computations 
for the (A) = case. We show that the theory has a ghost, and that the linearized perturbations diverge where the 
total mass of the vector vanishes. In Section |V] we study the case of nonvanishing vev. We first discuss the simplest 



Notice, however, that for a vector curvaton the vev cannot be exactly zero if one wants to realize a violation of statistical isotropy of 
the perturbations. Indeed, to have a violation of statistical isotropy, one needs to "single out" different direction(s); this is precisely 
provided by the vev of the vector. We follow the computations of 30] which are performed under the assumption of zero or negligible 
vev. For a nonvanishing vev, the model of [30ll belong to the second class of models just mentioned, and the ghost instability manifests 
itself in the second way just mentioned. 
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possibility of a vector and a cosmological constant. For this case, we both compute the quadratic action, and solve the 
linearized equations of the perturbations. Once again we verify that the modes diverge precisely when one eigenvalue 
of the kinetic matrix K crosses zero. We then study the case of one vector and one inflaton [2^, and the simplest 
realization of vector inflation [28| . For these models we compute the kinetic term of the physical perturbations, and 
we study how the eigenvalues evolve in time. Conclusions and discussions are given in Section IVTl 

II. REVIEW OF SOME MODELS WITH A RA^ TERM 

In this Section, we briefly review the reasons for introducing a nonminimal coupling of a vector field to the curvature. 
We start from the quadratic action of a vector field with a generic time-dependent mass: 



S= d 



74 „ 
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(2) 



leading the equations of motion 



^ ■.d^{y/^F'"') +M'^A'' =0 (3) 



Moreover, due to its antisymmetry, the field strength F^i, satisfies the identity 

dp Fp, + Fp,p + dp F,p = (4) 

Turner and Widrow [2^ suggested a mechanism for the generation of primordial magnetic fields during inflation, 
starting from the action ([2]). * In this mechanism, the vector field Ap is the electromagnetic field, and the mass 
term is proportional to the scalar curvature R. Following [27| . we use conformal time 77, defined by the line element 
ds^ = a? [rf) [—drf + da?^) in the present discussion, so that the electric and magnetic fields are (£123 = +1) 

Fio ^ , Ftj = e^fc Bk (5) 

The two equations ([3]) and (|4]) can be easily combined into an equation for the magnetic field (27| : 

(92 - 9| + M^) (^a^B^ = (6) 

From this equation, we find the well known result that, in the massless case (M^ = 0), and in the large wavelength 
limit (i.e. negligible spatial gradient) , the amplitude of a magnetic field decreases as oc a~^. Correspondingly, its 
energy density decreases as (x a"**. Consider instead the mass term 

= ^R = 6^^ (7) 

a-' 

where ^ is a constant, and prime denotes derivative wrt rj . The equation of motion for the magnetic field then becomes 

d-^,-dl + 6^^^{a'B)=0 (8) 

During inflation, a = — 1/ {H rf) where the Hubble rate H is nearly constant; inserting this into the equation of 
motion (and treating H as constant), we find that, in the large wavelength regime, the energy density of the magnetic 
field behaves as 

oc 52 oc a-s+vrac (9) 

Therefore, for ^ < - corresponding to a negative IvP in the action ([2]) - is less affected by the expansion with 
respect to the massless case. 



More precisely, ref. |27ll studied a more general action with a quadratic term of the type i?^i,Af'A'^ also included; we disregard this term 
in the present work. 



5 



It is instructive to compare the behavior of the magnetic field with respect to that of a massless scalar field coupled 
to the curvature, characterized by the action 



This leads to the equation of motion 



a^-9i+(66-l)- (a0) = O 



(10) 



(11) 



By comparing this equation with the analogous expression for the vector field ^, it has been noted [23| that a vector 
field with = —1/6 behaves analogously to a scalar field minimally coupled to the curvature {^s — 0). Conversely, 
the standard vector field ^ = is analogous to a conformally coupled scalar, £_s = 1/6 (therefore, no magnetic field is 
produced by the inflationary expansion in the standard case). 

This analogy has been recently exploited in [28| , that proposed a mechanism of inflation driven by a combination 
of N nonminimally coupled vector flelds. The mass term of the vectors comprises of the coupling to the curvature 
plus a constant term: 



N r 



5* = 



(12) 



The simplest case, and the one which has been most studied in |28| . is characterized by three mutually orthogonal 
vectors with equal vev: 



(AW) = (0, A 0, 0) , (42)) ^ (0, 0, A 0) , (43)) ^ (0, 0, 0, ^) , A = Mpa{t)B{t) 



(13) 



These sources allow for a FRW background, controlled by the equations of motion (we switch to physical time, and 
we denote by a dot derivative with respect to it) 



52 
2 



1 



1 + 6$ 



HB 2B+HB 



B + iHB + m^B = -{l + QClB[H + 2H^^ (14) 

For ^ — —1/6, and upon the identification B — (f>/ (-\/3Mp), we recover the same equations as those of chaotic inflation 
driven by a minimally coupled scalar field (j). 

Another compelling feature of the proposal of [2^ is that it can naturally give a small violation of isotropy. Indeed, 
for a large number N of vectors with random orientations and vev, one expects an almost isotropic expansion, with a 

deviation AH/H = O ^I/a/ZV^ between the expansion rates of the different directions. Ref. [29] provides a slightly 

different mechanism in which a single nonminimally couple vector breaks the isotropy in one spatial direction, while 
a scalar field with greater energy density is responsible for the overall nearly isotropic expansion. 

Vector fields with nonminimal coupling to the curvature have also been recently employed for the generation of a 
nearly scale invariant spectrum of perturbations [13] ■ Consider the action for a single field (N — 1) and with 
^ = 1/6. Following the discussion of .30 ], we compute the evolution of the perturbations 6Af^ assuming that the vev 
(A^) can be neglected (see the remark we made about this in the Introduction). In this way, the perturbations of the 
vector do not mix with those of the metric at the linearized level. Moreover, with a negligible vector vev, we can study 
the evolution of SA^^ in an unperturbed FRW background. In the present discussion, we only consider the transverse 

components of the perturbations, 6A^ = (0, SA^^ , with d, SAf = , since, as we will show in the following Section, 

there are serious problems with the longitudinal mode in these models. Due to the ^ = —1/6 choice, the equation for 
this mode is identical to that of a minimally coupled curvaton scalar field: 



sat^o 



(15) 



^ It is appropriate to compare the behavior of B with that of a<f>, since, due to the different structures of the kinetic terms, these are 
the canonically normahzed fields in the two cases. 
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We proceed as in the scalar curvaton case. For simplicity, we assume a dS background; in momentum space, the 
Fourier transform of 6Aj obeys to the equation 

{d', + k^-'-f)s^I-0 , M^^l-^ (16) 

Among the two solutions of this equation, we chose the one that reduces to the adiabatic vacuum at early times, when 
the mode is deeply inside the horizon: 

SAj^^VW\H^;i^i\k,\) (17) 

where H^^'' is the Hankel function of the first kind. Indeed, up to an irrelevant phase, this solution reduces to 
e~^'^^/V2k for \kr]\ ^ 1 . In the opposite late time / super horizon regime (\kri\ ^1), from the expansion of the 
Hankel function, we then find the power spectrum for the transverse modes [30 1 

oc fc3 \SAj\^ cx = (18) 

Namely m <^ H provides a small departure from scale invariance. We remark that this result follows from the 
solution (jl7p . Indeed, eq. (jl6p has two solutions with two undetermined (and, in principle, /c— dependent) integration 
constants. As it is customary, we chose the linear combinations of these solutions which reduces to the adiabatic 
vacuum in the early time sub-horizon regime. The phenonienological prediction (|18p crucially relies on the fact that 
the theory must be under control in this regime. 

For all the models we have reviewed, the discussion presented in this Section disregards the role of the longitudinal 
polarization of the vector field(s). In the next Sections we show that, for all the models discussed, this mode turns 
out to be a ghost. 



III. IDENTIFICATION OF GHOSTS, AND THEIR ASSOCIATED INSTABILITY 

In the next two Sections, we will see that the models described in the previous Section have ghosts. We outline here 
the method we employ to find the ghosts, and the instability associated with them. The ghosts are among the physical 
excitations of the background geometries of the various models. Therefore, we need to compute the spectrum of these 
theories. To do so, we perturb the background solutions of a given model, and we expand the action at quadratic 
order in the perturbations. This is the free action for the perturbations (interaction among these fields come from 
expanding the initial action to higher orders), and the spectrum follows from the diagonalization of this action. 

There are two main issues in this computation. The first one is associated to the gauge freedom that may be 
present in a theory. A gauge theory is in a sense a redundant formulation of a physical system, since different field 
configurations, related to each other by a (nonsingular) gauge transformation, describe the same physics. In the 
present context, the gauge freedom is the one associated with general coordinate transformations. We can expand a 
gauge transformation in a transformation acting on the background plus a transformation acting on the perturbations 
(loosely speaking, we decompose any gauge transformation into a "big transformation" , affecting the background, 
plus a "small" transformation, affecting the perturbations on a given background). In any explicit computation, one 
typically fixes the gauge freedom for the background and the perturbations separately. Concerning the latter step, one 
can choose a gauge that fixes the freedom completely (for instance, one can impose that some perturbations vanish; 
one needs to show that this choice can always be done, and that there is no residual freedom left); equivalently, 
one can find a set of gauge invariant linear combinations that do not change under the gauge transformation. Since 
the theory is gauge invariant, it must be possible to write down the equations of motion only in terms of the gauge 
invariant combinations, and solve those expressions. This is the method that we adopt in our explicit computations. 

The second issue is that, even after removing the gauge redundancy, the remaining perturbations do not all neces- 
sarily describe dynamical degrees of freedom (e.g., physically propagating excitations). Modes that do not correspond 
to dynamical excitations enter in the action without time derivatives (up to boundary terms, which can be disregarded 
in a theory without boundaries). Example of nondynamical degrees of freedom are for instance the SAq component 
of a vector, or the Sg^o components of the metric. In the standard case (i.e., for £ D F^^F^^" , and C D i?), they 
enter in the action without time derivatives. In general, the quadratic action for the perturbations around a given 
background is formally of the type (in momentum space) 



5*= / d'kdt 



^ij ^3 



(b,, N* i> + h.c.) + N* N, + (d., Y* Y, + h.c.) + e., Y* Y, + (/,, N* Y, + h.c.)] (19) 
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where Yi are the dynamical modes, and Ni the nondynamical ones. The coefficients , . . . , fij depend on baclcground 
quantities, and, for any given background solution, are functions of time. From the action (|19p . we find the linearized 
equations 



SS 
5Y* 

SS 
SN* 



= 
= 



Y, 



TV, =0 (20) 
(21) 



We see explicitly that the equations of motion for the nondynamical modes are algebraic in them (they are constraint 
equations). As a consequence, the nondynamical modes do not introduce additional degrees of freedom, but are 
completely determined by the dynamical ones. We can obtain a system of differential equations for only the dynamical 
modes by inserting the solutions of eqs. (PT|) into equations (PO)) : 



K Y 



where 



+ (Ay - h.c.) 



Y, 



A,: 



Y,=0 



(22) 



K., 



U — 
2 _ 
ij ~ 



y ^)km f"^i 



(23) 



To obtain instead the quadratic action for the dynamical modes, we insert the solutions of eqs. (j2ip into the action 



s 



dPkdt 



Y* K,, K,- 



'Y* 



A,,Yj 



+ h.c 



Y 



(24) 



The extremization of this action provides exactly the above equations (|22p . 

Eq. (pij) and ([^ - or, equivalently, eqs. and (pij) - are the equations of motion for the perturbations. We can 
also obtain these equations in a different but equivalent form by simply expanding at first order in the perturbations 
the general equations of motion of the model. In general, this last procedure is technically simpler than expanding the 
action at second order in the perturbations, and extremizing it. However, while an action provides the equations of 
motion, in general the equations of motion do not fully determine the action. The quadratic action of the perturbations 
is necessary to quantize the perturbations, and provide their initial conditions. 

The procedure just outlined formalizes the steps that are done in the standard case of scalar field inflation on a 
FRW background. Moreover, it extends it to an arbitrary number of dynamical and nondynamical fields. In the 
standard case, there are 5 scalar modes in the perturbations of the metric and of the scalar field. Two of them are 
eliminated by gauge fixing (or, equivalently, by the use of gauge invariant variables). Only 1 out of the 3 remaining 
modes is dynamical. While one can show that the action can be written solely in terms of gauge invariant quantities 
without using the equations of motion, the constraint equations are used to eliminate the nondynamical modes (for 
instance, the constraint equation (10.39) is used to obtain the final action (10.59) in [l^]). The final step in the 
computation is to normalize this final mode so that the kinetic term of this final action is canonical. 

For a general problem, with more than one dynamical mode in the system, one needs to diagonalize the kinetic 
matrix K. If one eigenvalue of K is negative, the corresponding eigenmode is a ghost (a field with negative energy). 
This signals an instability of the vacuum, which decays in ghost- nonghost excitations with a rate which diverges in the 
UV. A second type of instability takes place whenever, due to the cosmological background evolution, an eigenvalue 
of K crosses zero, and K is noninvertible. Denote by one of the moments at which this happens. The equations of 
motion become singular for t t^, , since the second derivative of the modes diverge in this limit. We expect that. 



In the standard case, the canonically normalized scalar perturbation is the Mukhanov-Sasaki variable v Issll . We explicitly verified 
that, in the standard case of scalar field inflation, the action H24| l obtained with the procedure described here coincides with the action 

(10.59) of m. 
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correspondingly, the solutions of the linearized system diverge for t ^ . In the next two Sections, we show that 
both types of instability occur in the models we are studying. 



IV. GHOST INSTABILITY FOR {A^) = 

We start our analysis from the case in which the vector field has no vev. This is the case for the computations 
presented in [23, [lOil to study the evolution of the modes. On a technical level, the assumption of zero vev drastically 
simplifies the computation. Indeed, the actions we are studying (cf. eqs. ^ and p^ ). are quadratic in the vectors; 
therefore, for (A^) = 0, the perturbations of the vector(s) are decoupled from those of any other field at the linearized 
level. Therefore, the action for the vector field already starts at quadratic order in the perturbations, and we can 
simply study the evolution of SA^ in an unperturbed FRW background. The more involved case of (A^) is studied 
in the following Section. 

Consider the action 



S 



1f^--m^ A^ 

4 2 



M 



-f2 



R 
6" 



(25) 



in the FRW background, ds^ ~ —dt^ + a [t) dx^ . Parametrize the fluctuations of the vector field as A^ ~ 
(ao, 9ia + af), where i — 1,2,3, and where the modes aj are transverse {diof — 0). These two polarizations 
are decoupled from the modes {ao, a} (and also from each other) at the linearized level, and have been already 
studied in the previous Section (starting from eq. (fT5|) ). As we saw, they are well behaved in all regimes. Here we 
study the nature and the evolution of the other two perturbations. 

It is well known that a massive vector has three physical degrees of freedom. The two transverse ones are encoded in 
aj . Therefore, the two perturbations ao and a encode only one physical mode, which is the longitudinal polarization 
of the massive vector. This can be immediately seen from the equations of motion for the two perturbations following 
from (P5|) . which, in Fourier space, read ^'^ 



ao 



p2 _^ f^2 

(3p2 + Af 2) H 



M2 



a = 



(26) 
(27) 



where p = fc/a is the physical momentum of the perturbation, and, as usual, H = a/ a. While eq. (j27p is a second 
order differential equation, eq. l|26|) is an algebraic equation in ao. Therefore, ao does not introduce additional degrees 
of freedom, but it is completely determined once a is known (compare these equations with the formal set of equations 
(|20| and ([2T|)). 

We can also see this from the quadratic action for the perturbations. Inserting the decomposition A^^ = (ao, di a) 
in (j25[) . and Fourier transforming the spatial coordinates, we find 



^~ 2 



dtd^k [p 



2| • |2 

'a 



- /Af>|2 - p2 a + c.) + (p2 + Af2) 



(28) 



The mode ao enters in the action without time derivatives, confirming that it is a nondynamical mode. We integrate 
it out: we compute its equation of motion from this action (namely, eq. ()27|) given above), we solve it, and we insert 



To integrate the system I I22I I. one needs to invert the matrix K\ in this way one obtains a system of equations of the form Yi = . . . , 
where the right hand sides contain only zero or first time derivatives of the unknowns. This system can be then integrated numerically. 
The inversion of K can be performed at all times t <tf, and the resulting equations are regular until t*; however, the expression for 
the second derivatives blows up for t — > t*, and we expect that the solutions diverge at this moment. This is precisely what happens in 
the two examples below in which we explicitly solve the linearized equations. This is also what happened in our previous study 31] of 
the model 26] which presents an identical instability. 

In the mechanism of |27l |. the full action in not quadratic in A^, since the photon is obviously coupled to the current of charged fields. 
This coupling is important from reheating on, when the conductivity of the charged particles become high, and affects the evolution of 
the magnetic field. However, the presence of these fields can be disregarded during the inflationary stage, as done in |27l |. 
The first equation is &S/&Aq ; the second equation has been obtained by combining 55/(5 Aq and SS/&Ai . 
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the solution back into (1281) . This leads to the action for the longitudinal mode: 
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S=ljdtd^kp^M^[^-^^-\an (29) 



The extremization of this action reproduces the equation of motion ^7\ . We remark that the lagrangian in 
is proportional to the term, and that, for > |M^|, the longitudinal mode is a ghost (i.e. a field with negative 



energy, not simply a tachyon) whenever < |25l |3l| . In our previous work 3l| we provided two additional proofs 
of the fact that < leads to a ghost. The first is based on the direct computation of the propagator (we find 
that one residue at the pole of the propagator has the "wrong" sign for AP < 0); the second on the Stuckelberg 
formalism (which is a convenient way to study the different polarizations of a massive vector or graviton; we rewrite 
= Aj^ + ^fJ,(j)/^/\M^\, with d^" = 0; it is then immediate to see that the field is a ghost for < ). 
The presence of a ghost signals the instability of the vacuum of the model, due to its (UV- divergent) decay in 
ghost-nonghost excitations. In addition, the equation of motion (j27p indicates that the system may be unstable 
already at the linearized level, whenever = + M"^ or vanish. Let us first discuss when this occurs. During 
inflation, ^ R for the models we are studying. As a consequence, during inflation, AP — m? — i?/6 ~ —2H^ < 0, 
while w'^ = p'^ + A-P ~ — 2H^ goes from positive to negative. We denote by t^i the moment at which lo^ vanishes; 
we note that this happens when the mode is close to horizon crossing, and that lo^ remains negative from t^^i until a 
time well after the end of inflation (since the mode is well outside the horizon, p H, when inflation ends). After 
inflation, R/6 — O (H^) decreases as the universe expands, and it eventually drops below m^. At this moment, which 
we denote by tM, the mass term goes from negative to positive. Therefore, = + AP > for all times 
t > tM ■ Since we saw that < at the end of inflation, there is a moment after the end of inflation, and before Im, 
at which vanishes for the second time. We denote this moment by ti^2- In summary, vanishes ai t = Im, while 
Lu'^ vanishes a.t t — t^ji, ^^2- Denoting by icnd the moment at which inflation finishes, and by Iq the present time, we 
have tail < iond < iw2 < tn < to ■ 

To see whether the linearized system diverges when either or A/P vanish, we study the equation of motion (j27p 
for a finite interval of time close to t^ji (i = 1,2) or t^i. We assume that the equation of state w of the source driving 
the expansion can be treated as constant in this interval (which is certainly true, provided the interval is not too 
extended). We also assume that —1 < w < 1/3. This is a rather general assumption, since it includes the following 
cases: inflation (with w > —1), coherent oscillations of the inflaton after inflation (which, for a quadratic inflaton 
potential, give an average u; = 0), matter domination (if is sufficiently small, so that AP vanishes at this stage), 
and also radiation domination (for which the equation of state is slightly smaller than 1/3, due to the masses of the 
particles in the thermal bath, or the thermal trace anomaly). 

Behavior of the linearized system for p^ + AP . The scale factor evolves as 



I \ 3(1 + ™) 

a = a^A—] (30) 



were we recall that t^i is either of the times at which lo^ — 0, and a^^ is the value of the scale factor at that time. 
The mass and frequency squared are given by 

,,r2 2 R 2 2 1 - 3w 

M ~ m — — = m ——7 ^ 31 

6 9^2 (l + u;)2 ^ > 



. -,^ + M^.,;.|^) +..."-5-^ (32) 



2 l-iw 

\2 



We remark that, in going from the action ||28|I to the action (|29| l. we are precisely following the general procedure that we outlined at 
a formal level in the previous Section. 

This equation, in conformal time, was already given in |30| : this confirms our algebra. The action for the longitudinal mode could be 
also written starting from this equation, up to an overall factor. The procedure described here provides the complete action. We also 
note that the Ref. 30] did not solved the equation at the moment in which the total mass of the vector vanishes; as we show here, the 
solution diverges at this moment. 

One may imagine that m is so small, so that tm has not occurred yet. In this case the longitudinal polarization is still a ghost today; 
we disregard this possibility, due to the stringent limits on theories with ghosts found by |i35j and discussed in the Introduction. 
The thermal trace anomaly is relevant for temperatures greater than the QCD phase transition, and gives 1/3 — «) = O (10~^) [39|l . 
Mass thresholds give 1/3 — w = O (lO~^) when the temperature is close to the mass of a particle [4Cl| . For temperatures below any 
Standard Model particle, the departure from jd = 1/3, and _R = can be neglected. However, as R drops to zero, there will be a moment 
in which vanishes; we denote by w the value of the equation of state at this moment. 
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where p^i is the value of the physical momentum at t^i. We then find 

%/2 VI - 3w 



w (Ci) = 



3 (1 + w') 



(33) 



We insert these expression for the momentum and the total mass in eq. (|27|). and we Taylor expand the resulting 
expression for t ~ t^ii- We find: 



t-t. 



+ C (t-t^,) awO 



where 



Eq. dSl]) is solved by 



^=V^T-^ [3m^(l + ^)+pL(l + 3u;)] 



^1/3 



(34) 



(35) 



(36) 



where Ci^2 are two integration constants, and Ai' and Bi' the derivatives of the Airy functions Ai and Bi , respectively. 



These two solutions are regular at t = ti^i, where they have the expansion series Ai' , Bi' 
the linearized term is absent, we find that oc a/ (t — ti^i) also remains finite as t = ti^ 

Behavior of the linearized system for . The scale factor evolves as 



const. + (< — 



Since 



a = aM 



t \ 3(1 + ™) 



were we recall that Im is the times at which — 0, given by 

V2 VI - 3w 
tM = -5 — 

6171 1 + W 

(cf. eq. (PT|)). and cm is the value of the scale factor at tM- We Taylor expand eq. l|7f)) for t ~ tjy/: 



a - 







t-tM 

where pM is the value of the physical momentum of the mode at tM ■ This equation is integrated to give 

a « Ci Jo (pm (tM - t)) + C2 Yq {pm (tM - t)) 



(37) 



(38) 



(39) 



(40) 



where Jq, Yq are, respectively, the Bessel functions of the first and second kinds of order 0, and Ci,2 are integration 
constants. While the Jq solution is regular at i = ^m, the Yq solution has a logarithmic divergence. Correspondingly, 
the mode ao exhibits a linear divergence, as can be seen from equation ([26|) . 

The only way to avoid the singularity is to arrange the initial conditions so that C'2 = 0. This must be done for 
every mode (namely, for any comoving momentum k) and for both the real and imaginary parts of the perturbations. 
We regard this as a completely unnatural assumption, since there is no reason why the initial conditions (set during 
infiation) should "know" about the singularity which is to occur later on when vanishes. 

In conclusion, the solutions of the linearized system remain finite when = + — 0, while they diverge 



when AP = 0. To verify this behavior, we performed a numerical evolution of eq. ((27)) for the specific case of w = , 
Pm = 2 . We then have tt^2 — 0.058/m and tM — 0.47/m . We set the initial conditions a = 1, d = at iAf/lOO . 
We plot in Figure [T] the resulting evolution of a and ao . We see that the system is regular at t^2 , while it diverges at 
tM (we verified that ao diverges linearly, and, correspondingly, a diverges logarithmically). 

It is worth noting that the solutions diverge when the kinetic term of the longitudinal mode vanishes, cf. eq. (|29p . 
On the basis of what we wrote at the end of the previous Section, one should expect that this happens also for the 
cases in which the vector field(s) has a nonvanishing expectation value. We verify this explicitly in Section fVBl for 
the simplest possible case of a vector field and a cosmological constant. 



The parameters of the evolution have no particular relevance, and have been chosen only for illustrative purposes. We have verified that 
other choices of parameters also confirm the behavior we have obtained analytically. 
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FIG. 1: Behavior of the hnearized solutions representing the longitudinal mode. The system remains finite at the point tu:i, 
when the frequency vanishes, but it diverges at Im, when the mass term of the vector vanishes. See the main text for the 
parameters chosen. 

V. GHOST INSTABILITY FOR (A,,) / 

As we mentioned, when the vector field(s) has a vev, its perturbations mix with those of the metric already at 
the linearized level. The computation then becomes significantly harder than that discussed in the previous Section 
for the case of vanishing vev. Besides the increased number of perturbations, the spatial vev of the vector (s) breaks 
the isotropy of the background, so that, in general, one needs to study the cosmological perturbations of a non FRW 
space. 

We start by studying the case of a single vector field with a vev along one spatial direction. The background has 
a 2d isotropy in the other two directions. In Subsection IV Al we review the formalism for dealing with the coupled 
system of perturbations of the vector and the metric in this context (l^ . [sij . By exploiting the symmetries of the 
background, the perturbations can be classified in two distinct subsets, which are decoupled from each other at the 
linearized level. We verified that one subset of perturbations, which includes the transverse vector modes, is well 
behaved, and, for brevity, we disregard it in the following discussions. The other subset, which we discuss in details, 
includes the longitudinal polarization of the vector. 

In Subsection I V Bl we apply this formalism to the simplest case of a single vector field plus a cosmological constant. 
In Subsection IV CI wc study the case in which the cosmological constant is replaced by a scalar inflaton, which is the 
model proposed in |2^. Finally, in Subsection IV Dl we study the model of vector inflation [2^ (for which the formalism 
of Subsection IV Al does not apply; for this reason, this Subsection is a self-contained study). 

A. General formalism for a vector field in a 2d isotropic background 

We review the general formalism for the study of perturbations [l^, [3l| for the case of a vector field with a 
nonvanishing vev along one spatial direction: 




(A^) - (0, Mpa{t) Bi(t), 0, 0) 
ds^ = -df + [t) dx^ + I? [t) [dy^ + dz^) 



(41) 
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The background geometry is a Bianchi-I background with 2d isotropy. The complete set of metric and vector field 
perturbations can be written in a way which exploits the residual 2d-isotropy of the system: 



-2$ adix b{d,B + B,) 

-2a^^ abdi {diB + B,^ 

9 (-2S (5y -2didjE- o/e. 



(ao, ai, dia + ai) 



(42) 



where i,j = 2,3 span the isotropic 2d-plane. The perturbations {<!>, x, B, B, S, E, cn, ct} are 2d scalar 
perturbations, comprising one degree of freedom (d.o.f) each, and Bi, Ei, a^} are 2d vector modes which satisfy 
the transversality condition (5^ Bj = • • • = 0). Therefore, the 2d vector perturbations also comprise one d.o.f. each. 
Contrary to the 3d decomposition, there are no tensor modes, since there are no degrees of freedom left in a symmetric 
2x2 transverse and traceless matrix. The vector and scalar modes are decoupled at the linearized level (namely, in 
the linearized system of equations for the perturbations, or, equivalently, in their action, obtained by expanding the 
action for the model at quadratic order in the perturbations). 

We need to eliminate the redundancy associated with the freedom of general coordinate transformations. We can 
do so by either choosing a gauge that completely fixes this freedom, as done in [l^ , or by combining the perturbations 
in gauge invariant modes [31|]. We choose this second method, since it provides a nontrivial check on our algebra 
(since one needs to show that all the equations, and the action, can be written solely in terms of these combinations) . 
Among different equivalent choices, we use the gauge invariant combinations (sij 



B 



Mr, 



+ (b+^-e 

Hb a \ a 



B- 



1 



P «2 



1 



Y. + bE 



ai 



S-a ( - (b + -E 

Bi + Ha Bi 



ai+a Mp " ' " E - 6 Mp Bi d'f [B 



-la 

a 



Hb 



a~bMpBidi I B + -E 



E 



do = -9i 



aa-aMpBidi ( - [B + -E 



(43) 



(where the dot and buhet denote time derivative, Ha ^ a/a, Hi, ~ b/b , and 9|n = 9| + d^) in the 2d scalar sector, and 

Hi = Bi -f- b Ei 



B, = a{ B, 



■E, 



Mr, 



(44) 



in the 2d vector sector. If there is also one scalar field. 



there is the additional 2d gauge invariant scalar mode 



Hb 



(45) 



We note that the gauge invariant modes have the following mass dimensions: 



B, 



-1. 



B, 



0, 



1, 



B 



[x] = [«o] = 2. 

The above gauge invariant combinations of the metric perturbations do not reduce to the ones which are commonly 
given in the literature fSTt] in the limit of isotropic background. However, we explicitly verified in [3lj that our choice 
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and the more conventional one are equivalent. Our choice is motivated by the fact that it immediately identifies the 
nondynamical modes of the system: one can choose the gauge B = E = T, = Ei = 0, in which the gauge invariant 
combinations Sqq^ and SAq reduce to the corresponding non dynamical perturbations ^^q^ and Aq. These modes 
enter in the action without time derivatives, due to the specific form of the "kinetic terms" and R. Therefore, 
we expect that also the gauge invariant combinations SgQ^ and SAq are nondynamical, as the explicit computations 
reported in the next Subsections confirm. 

As we mentioned, the two subsets of modes (2d scalars vs 2d vectors) are decoupled from each other at the linearized 
level, and can be studied separately. We verified that, for the models of our interest, the 2d vector modes do not 
develop any instability. For brevity, we do not report these computations here, and we focus on the more problematic 
2d scalar modes. In the next two Subsections we study the evolution of this system for the case of a single vector plus 
a cosmological constant or a scalar inflaton. 



B. One vector plus a cosmological constant 



We study the simplest model in which a nonminimally coupled vector field with a spatial vev provides an anisotropic 
expansion. Besides the vector field, there is a vacuum energy Vq which is responsible for an overall accelerated 
expansion. The Subsection is divided into three parts. We first present the model, and discuss the background 
evolution. We then solve the linearized system of equations for the perturbations, and find that it diverges at some 
point close to horizon crossing. We finally compute the kinetic term K of the quadratic action of the perturbations, 
and show (i) that the divergence of the linearized system takes place precisely when one eigenvalue of K vanishes, 
and (ii) that one of the perturbations is a ghost for some moment of time. 



1. The model and the background solution 



The action of the model 



S 



4 ^ " 2 I 6 



(46) 



gives the equations of motion 



Gn 1/ — 



IvP 



1 



m^^^ \ A^A^^ ^m' A^ A° g^,^ 



1 



Rt.v - 2 ^.9m- ) ^" ^" - « (fff- n - V,) A^ 



:9u 



6 

A"^ = 



(47) 



If we assume that the background solutions depend only on time, then the v component of the last of (j47p gives 
A^ — Q . For a homogeneous background, the vev of the vector is everywhere in the same direction; we choose the 
coordinates such that this direction coincides with the x— axis. We look for background solutions of the form (I4ip . 
and we define 



a 
a 



H, 



rr Ha + 2 Hi, Hf, — Ha 
H = , h = 



(48) 



This was the choice made in [12| |. and it fixes completely the freedom associated to general coordinate transformations. 
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Namely, H is the average expansion rate, while h parametrizes the departure from isotropy. The regime of small 
isotropy corresponds to /i <C i? . In terms of these quantities, and for the background (|¥T|) . equations (|T7|) give 



h + 3Hh 



Vo 


1 


3 Af 2 


+ 6 




H 





-hBiBi + -Bf (m^ - mh + 5/1^) 
3 6 

+ ^Bl + i (2i? - bh) BiBi + i Uh^ - ^Hh + bh^ - ) 



2ij + 3i/2 + 3/1^ ^ - - ]-Bi \Bx + {3H - 2h) + \ {iHh - + m^) Bl 



Bi+3H Bi + {m^ -5h^ -2hH ~2h] Bi = 



(49) 

One of the last three equations can be obtained from the other equations in (|49p due to a nontrivial Bianchi identity. 
We see from the second of (j49|) that the anisotropy of the background is indeed supported by the vector vev Bi. 

Below, when we study the perturbations of this model, we use the exact background equations (j49p . However, only 
for this discussion, we present the inflationary solution of (|49p in the slow roll approximation (slow motion of Bi 
towards zero, so that the anisotropy is not quickly damped away) and for small anisotropy. Specifically, we disregard 
the time derivatives of H, h, and look for expansion series solutions of the type 



H^Ho+chBI + {Bf) 



ChB't + 0(Bf) , Bi~cbHoBi , with Bi, cb < 1 



(50) 



We insert these expression in the last equation of and expand in Bi. This gives 3cbH^ + c\ Hq + — . 

Since, cb <C 1, this gives cb — —-ttti- Therefore, slow roll requires m ^ Ho . We then insert (l50l) in the first of (1491). 

expand in and neglect m vs. Hq. We find Hq = \/VqI {yfZM^ , cu = m^/ {\2Hq) . Finally, we insert ([50)1 in 
the second and third of (|49p . expand in _Bi, and again neglect m vs. H^. These equations then become identical, and 
give Ch = Hq/2>. Therefore, 

H = H, + -^^Bl + 0{Br) , h=^Bl + 0{Bl) , B, = -^B, + 0{Bl) (51) 

where Hq is the expansion rate in absence of the vector field. 

Before moving to the study of the perturbation, we note that the mass parameter m of the vector field must satisfy 



Bl Hi < < Hi 



(52) 



The second condition is due to the slow roll requirement, as we already pointed out. The first condition comes from 
requiring that the "mass term" for Bi in the last of (I49p is positive, so that Bi evolves towards the origin (so that 
the anisotropy decreases during inflation rather than increasing). 



2. Instability from the linearized equations 



We now expand the fleld equations (|47p at first order in the 2d scalar perturbations of the metric and the vector 
field. 21 We denote the resulting linearized equations as 



We then Fourier transform these equations. Namely for each mode, we have 

5{x) = [ dPk5{k) e-*fci^-*'=^2«-^fcT3^ 

(27r) J 







(53) 



(54) 



The first of I I49I I is the 00 Einstein equation; the second is the hnear combination of the (ll)-(22) Einstein equations; the third equation 
is the hnear combination of the (11) + 2 X (22) Einstein equations . Finally, the fourth equation is the x— component of the equations 
for the vector field (the 33 Einstein equation coincides with the 22 one, while the remaining equations are trivial). 

As we mentioned, the 2d vector perturbations are decoupled from the 2d scalar ones; we have verified that the 2d vectors do not develop 
any instability; for brevity, we do not present this computations here. 
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where S denotes any of the perturbations (we denote the perturbations in coordinate and momentum space with 
the same symbol). Modes with different comoving momenta are decoupled from each other at the linearized level. 

We denote by /cl the component of the comoving momentum in the x-direction, and by fc t the component of the 
momentum in the y — z plane. The full comoving momentum is then given hy k'^ ~ kf^ + k . The physical momentum 
is instead 



2 2 I — >2 

P ^Pl + Pt 



k 'j-i 



(55) 



The explicit expressions of the linearized equations for the scalar sector are given in equations (jAl|l of Appendix [Al 
We could express these equations solely in terms of the gauge invariant modes defined in (j43p . This is a nontrivial 
check on our algebra. Here, we disregard some of the equations in (|A1[) which can be obtained from the remaining 
ones (due to Bianchi identities). The (complete) set of independent linearized equations which we choose to integrate 
is 



Eqo 



Eq 



-01 



Eqoi 



Eqn 



Eq 



11 



1 



5? 



H 
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-BiBi 
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Pt 



2 Pt + 2pi 



6 



1 



Bt + 2 [BiH + Bi] Bih 



6 



~AhH ] Bt-4:hBiBi+Bi' 
X 



iH + h) + lBiBi 
6 



6 



1 fp^ 2 

--rn 



6 ^ 



B- 



6 



+ 4:hH] Bi + h B 



Bi 



-H -h] Bi + -Bi 
2 J ^ 2 

di = 



1 



- ( ai - * ) Bi + — [BiH -2Bih-2Bi] * 



6 



Bl 



Bl 



H + h 



Bl Bl 



6 + B( 



9 

Pt 



Pt 



-BiH + 2hBi + Bl ai 



6 + i3f 6 Bl , 

^ B Y '^•^aaa "0 - U 



- 1 - * + - Si di 



6 



^iHBi-2i3i]-['- + ^]h 



Bi 



12 



1 



7 



2 • 



H 



Bl 



h-^\Bi-^ 
B = 



1 Bl\ u 

2 + Tf 



1 

l + ^)H+-BiBi 



-H - - h] Bi + -Bi 
3 6/3 

*-7 fl + ^l X 



ai 



ai 



Bl 
6 



2hBi- H Bl - Bl] + {H -2h) ai 
1 



Pt_ 

pI 



Pt_ 

pI 



{H-2h)a+[HBi-2hBi + Bi) ^ - —V^^x 



^ „2 q '-^aoao 

Pl 



do = 



(di - 2 Bi^^ - (3 H Bl + 2 Bi^ Bi * + ^ (^-HBi+8hBi - Bi^ di + {M^^, - Pr Bf) * 



9 

— m 

4 

1 



Pt 



Bl di 



(3-Bf) H+\biBi 



Pt + 2pi _ 15 ^ 2 



---Bf B + -Btx+{2H-h) Bix 



2hBi- HBi- Bl] do = 



3/ii7 ) Bf + 3hBiBi -3 
9 1 , 



Pt 
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^h" --H^ + -Bl 
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Eqi 



Ai- r^Biii + 3HAi+ [^Bih-'^BiH - Bi) ^ + {Mc.,c., + Pt) 



Si 



Bi 



Bi-2Bihj $ + -3^ [B + x]+ao + ^ [p' - 6m^) $ + - {7 Bi h + Bi H - 3Bi] B+- {2H-h) Bix 
-p'^a + 2 {h + H) ao = 



Eq, : a + 3 {H -2h) A + {Maa 



pI) 



a-p\ai 



9 

Pi 

Pt 



BiH - 2Bih + Bi) B + ao + {2H - Ah) ao = 



(56) 



where T)^^, Vagaoi -A^**, ■M^an -A^qiqi, Maa depend on the background quantities and are expHcitly given in 
equations (jA2p and (|A3|) in Appendix [Xj In writing (j56p . we have also made use of the physical momenta defined in 
eq. ([551). 

We now solve the system of equations (155)) numerically. We start by noting that the first four equations in ([55)) 
contain at most a single time derivative for the perturbations, and do not contain any time derivative of the modes 
{$, X, B, do}; these are the nondynamical modes of the system, and these first four equations can also be obtained by 
extremizing the quadratic action for the perturbations (eq. ()63p below) with respect to these modes. These equations 
are precisely of the type (|2ip given above where we discussed the linearized equations at a formal level. Compare also 
with eq. ([^5)) for the nondynamical mode ao in the simpler case studied in Section ITV) 

One can choose to solve these four constraint equations for {$, x, B, do}, and insert the solutions in the remaining 
equations in ([56|) . In this way, we obtain a closed system of equations for the dynamical modes. A different but 
equivalent way to integrate the system ([55)) is to diiTerentiate the constraint equations, so as to obtain differential 
equations for the nondynamical modes too. Combined with the remaining equations in ()56|) . one then obtains a closed 
system that can be numerically integrated (the constraint equations need to be imposed as initial condition). In the 
first method mentioned, one obtains a system with fewer but more complicated equations. On the other hand, the 
second method results in a larger set of less complicated equations. Here, we adopt an "intermediate" method, which 
we have found to be convenient for the numerical integration. We solve the second of (|56p for x: 



Pi 
V 



(6 + Bf) B + 8B1 



10 <^o 



[{H -2h) Bi 



i-Bi] ai-4Bi 



[{H-2h) 



-4 (6 + B^) W !> - 4 Bi di + ABl * 



Bi ~2Bi] * 



(57) 



where the time dependent coefficients T>aoa(n ^ ^^id H are given in equations ()A2[) and ()A4[) of Appendix \X\ 

It is important to verify that the second of ([55)) can indeed be solved in terms of x', namely, that V ^ 0. The easiest 
way to verify this is to use the slow-roll solutions ()5ip in the expression (IA4[) for D, since we are working in a regime 
in which these slow-roll solutions are highly accurate. This gives 



V^6p^+ (p| - 24 + Urn^) Bf + {Bf] 



(58) 



We can integrate out x provided that V ^ . We see that V is positive at least as long as pt ^ 2 H Bi (when the two 
terms become equal, the negative term dominates the O (-Bi) expression, an the expression cancels with the positive 
O (Bi) term. At this point, the expansion series ([55)) breaks down). As we show below, the computation we are 
performing shows that the linearized perturbations blow up when p^ ~ \/V2 H / Bi . Since integrating out x is a step 
of this computation, our result can be trusted as long as pT > 2H Bi when p^ ~ \/V2 H / Bi. Namely, we consider 
modes for which px > (S^/a/S) pl at that moment. This is not a strong restriction, since Bi <C 1 . 

We insert the solution ()57p into the other equations in ()56p . and we differentiate the three remaining constraint 
equations (we stress that we do not lose any information, provided that these equations are imposed as initial 
conditions). In this way, we obtain a system of 6 differential equations, which can be formally written as 



di 
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B 

V do y 
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K44 '«45 't46 
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(59) 



V '^61 '^62 '^63 '^64 '^65 '^66 / 
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which correspond, respectively, to Eq^;^, Eq;^, Eq^, and to the dime derivatives of Eqop, Eqg^, and Eqp. 

The elements of the matrix depend on the background quantities. The functions on the right hand side of 

([55)1 fi are expressed as linear combinations of {^, ^, ai, di, a, a, B, ao} whose coefficients also depend on the 
background quantities. The explicit expressions for the matrix Ai^. and the coefficients fi are given in equations (|A8[) 
and (|A11|) of Appendix El 

We invert A^^i and integrate the system numerically. The determinant of A^^ vanishes at some given time. As we 
approach that time, diverges. Correspondingly, the numerical solutions of the system also diverge. 

By inserting the explicit expressions for the entries of given in Appendix [X[ we find: 

detA^, = g^ellf (/i + )' [72 -6Bf{3- 2I?„„„„ ) - (3 - 8P„„„„ )] (60) 



Using the expressions for 'Dagao and T) given in (jA2|) and (jA4p . we then find that the determinant vanishes when 
1 



2 _ 2 _ 

Pl — Pl* — u2 la , Tj2\ 



Bl {e + Bi 



18 (2m^ - AH^ +Ah^ + Bfj - (Gm^ - UH^ -72hH + 102 + IB^j Bf 
4H^ -28hH + 31 h^) Bf + 48 h Bi Bi + 4 {7h - 2H) Bl Bi I (61) 



This expression for p^^ is exact. Using the slow-roll solutions (|51[) . we obtain the approximate expression: 

, 6(2iJ^-TO2) 

pi, ^ '- + O (62) 

To confirm that the linearized modes indeed blow up at pl*, we integrate the system (|59p for a specific choice of 
parameters (the parameters we use have no particular relevance, and are chosen only for illustrative purposes). We 
choose a small anisotropy, {h/H).^^ — O {B\^^ ~ 10~^ , and an initial value for the momentum, ~ lO'^ifin, so 
that the mode is initially deeply inside the horizon. As we verify in Appendix |B1 the frequency for the modes is 
initially adiabatically evolving, so that we can choose the initial conditions for the dynamical modes 5*, di, d and 
their derivatives according to the adiabatic vacuum prescription. We then use the constraint equations (the first, the 
third, and the fourth of (|56p ). to provide the initial conditions for the three nondynamical modes $, _B, do of the 
system. The resulting evolution of the "relativistic Newtonian potential" $, and of two other 2d scalar modes are 
shown in Figure O We see that the modes indeed diverge at some given time. We verified that pl = Pl* at this 
moment. This is confirmed by the time evolution of one of the eigenvalue of the kinetic matrix of the dynamical 
perturbations (computed in lVB 3|) shown in the left panel of the Figure. We see that the linearized solutions blow up 
precisely when the system (|59p . or equivalently, the kinetic matrix, becomes singular. 



3. Ghost from the quadratic action 



In order to understand the physical reasons behind the instability we have just found, we compute the quadratic 
action for the perturbations. We expand the action ((^S)) at quadratic order in the perturbations l|42p. The resulting 
action can be expressed solely in terms of the gauge invariant modes defined in (I43|) - (|44p . This provides a nontrivial 
check on our algebra. We also find that the action is the sum of two separate parts, one involving only the 2d scalar 
modes, and one the 2d vectors. We disregard this second piece in the following discussion. The action for the 2d 
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FIG. 2: Results from a numerical simulation with — 0.1 //o -Bin ~ 0.1, PL,in = 400i/in, PT,in = 900//in, (corresponding 
to H/p ~ 10~^; the modes are initially in the adiabatic vacuum; only the final part of the evolution is shown in the two 
figures). Since H/p grows during infiation, we use this quantity as "time variable" in the Figure. Left panel: gauge invariant 
relativistic gravitational potential <1>. We show the real part in units of 'i'in. We also show the eigenvalue Ai of the kinetic 
matrix (multiplied by 3 x 10"" , so that it is visible in the figure) . We see that "1> diverges when Ai = . Right panel: real parts 
of the modes ao (in units of 100 Ho '^in) and a (in units of ^'in) . Also these modes (as all the modes of the system) diverge 
when Ai = . 
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BiH -2Bih + Bi] I do + h.c. ) - 
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(63) 



where the coefficients {I?>t*, I'qq, I'aiai, T^xx^ T^aaaoi ^xaol depend on tlie baclcground quantities (and, 

hence, are functions of time), and their expHcit forms are given in equations (|A2p of Appendix [Xl As a further 
check on our algebra, we explicitly verified that the extremization of this action with respect to the fields included in 
it reproduces the system of linearized equations 

In Section IIIII we outlined at a formal level the computation that we are now performing for the model (I46p . 
We expressed the quadratic action in eq. (fT9|) . where we distinguished between the dynamical fields Yi and the 
nondynamical ones Ni. The action is the explicit form of the action for the model we are studying. 
The variables x, B, do} are the nondynamical modes {Ni}, since they enter into the action ([M]) without time 
derivatives. Starting from the action (|63p . we integrate out the nondynamical modes, and compute the action for the 
dynamical modes, following the same steps that lead from eq. to eq. In practice, we read the coefficients 

Oil", . . . : fij from the action (j63p . by comparing it with the formal expression pop : we then compute the combinations 



K 



f^c ^ f (cf. eqs. (US])) that characterize the action of the dynamical 



A,; 



are 



modes. This computation is a straightforward algebraic manipulation; the resulting coefficients K^j, 
extremely lengthy, and we do not report them here. 

The instability emerged from the linearized equations can be understood by studying the kinetic matrix K . In Figure 
[3] we show the three eigenvalues of this matrix, for the same numerical evolution (i.e., for the same parameters and 
initial conditions) as the one leading to Figure O We see that the two eigenvalues A2,3 are always positive, indicating 
that the two corresponding eigenmodes are well behaved positive-energy fields. On the contrary, the eigenvalue Ai 
vanishes close to horizon crossing. The system of linearized equations becomes singular at this point (cf. the formal 
equations (|22p: they are singular if the matrix K is noninvertible), and the linearized solutions diverge. We also see 
from the Figure that the eigenvalue Ai is negative for some time after this moment. The corresponding eigenmode is 
a ghost in this time interval. 

Although the exact expression for K is rather lengthy, it is possible to obtain a simple expansion series for its 
determinant, in the sub- horizon / early time limit: 



det 



K 

a 6^ 



96 p6 



16p2p6 



O 



El 

p4 



(64) 



To obtain this expression, we first expanded the exact expression for the determinant in a power series of the momenta; 
we then simplified each term in the series by using the slow roll solutions (|51[) . and finally kept for each term only 
the leading expression for <C 1 (which corresponds to small anisotropy, since h/H = O {B\^ ). The terms of 
O {Hq/p'^^ are parametrically suppressed with respect to the second term in (|64p for Ho <C p. Therefore, the first two 
terms in (|64p provide an accurate approximation of the determinant in the whole sub-horizon regime. 

Eq. (|64p shows that the determinant is positive at sufficiently early times, and it then becomes negative in the 
later part of the sub-horizon regime. This confirms the behavior of Ai seen in Figure [3] (since A2,3 > 0, the sign 
of the determinant coincides with that of Ai). The determinant vanishes when the two leading terms in (j64p are 



(approximately) equal; this happens for ~ 6 



{2Hl 



value (|62p of pl at which the linearized system (j59p becomes singular. 



/ Bi . Not surprisingly, this is precisely the approximate 
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FIG. 3: Evolution of the eigenvalues of the kinetic matrix. The parameters and initial conditions are as in Figure [2] Due to 
the wide range spanned by the eigenvalues, in the y axis, we have used a linear scale inside the interval [—0.01, 0.01], and a 
logarithmic scale outside. 




FIG. 4: Determinant of the kinetic matrix, for the same choice of parameters and initial conditions as in the previous Figure. 
Compared with the previous figure, we show a close up of early times, around the point where the determinant vanishes. The 
black dashed curve shows the exact determinant, while the red curve shows the approximate expression given in eq. (|64|) . 



In Figure [H we compare the approximate expression for the determinant - the first two terms in (j64p - with the 
exact one (not reported here). We see that the approximated expression is extremely accurate in the range of our 
interest. 



C. One vector plus a scalar inflaton 



The model proposed in ref. f29| is very similar to the one we have studied in the previous Subsection, with the 
only difference that the vacuum energy that we have included is replaced there by a slowly rolling inflaton field. We 
expect that the same instability that we have found above is present also for this model. The study presented in 
this Subsection confirms this. We first briefly review the model and its slow roll solution. We then show that one of 
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the 2d scalar modes is a ghost. More precisely, as in the case studied above, the mode is well behaved in the deep 
UV regime, but it becomes a ghost close to horizon crossing. When this happens, the system of linearized equations 
becomes singular. 



1. The model and the background solution 



The model of ref. [29| is characterized by the action: 



S= d 



1 



1 



(65) 



The inflaton field cf) replaces the vacuum energy Vb that we considered in the action (|46|) . In this way, one can have 
a graceful exit from inflation. The field equations following from this action are 
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(66) 



where T^J^'' is defined as in (|T7)l . The background solution is again of the form (HH), giving 
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AP 

Bi +3HBi + (^m^ - dh^ - 2h H - 2h] Bi = 
4> + 3H^ + V' (0) =0 



(67) 



where we have used the same combinations of Einstein equations we have written in (|49p . 

We assume that the inflaton field is in the slow roll regime, 4> « —V'{(j))/3H , leading to inflation. For what concerns 
the evolution of the vector field, and of the anisotropy, we all the considerations done for the model (|46)) are valid 
also in the present case, with the only difference that the vacuum energy Vq is now replaced by the slowly decreasing 
potential energy of the inflaton. This leads to the slow roll and small anisotropy {Bi ^ 1) background evolution 



^^Bf + 0{Bt) , h^^Bf + 0{Bt) , B,^-^^B, + 0{Bf) , 



5Ho 



O (B!) 



Ho 



(68) 



2. Ghost instability from the quadratic action 

We now study the perturbations of the background solution. We again decompose the perturbations as in (I42p . 
disregard the well behaved system of 2d vector modes, and construct the gauge invariant combinations (j43p and (I45p . 
We proceed as in IV B 31 by expanding the action of the model at quadratic order in the 2d scalar modes, and by 
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rewriting it solely in terms of the gauge invariant modes (which provides a nontrivial check on our algebra). We find 
S2ds ^ \ I d^kdtab^ £2dS 
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where we have defined 
1 



3 I + 12ft2 - 12i?2 + + 3^ j -B? - 24ft Bf Bi + 2 {2H ~ 7ft) Bf Bi 



pI {18 + 3Bl + 2Bt) I ' A/2 
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In the action (j69p . we included only the terms that contribute to the kinetic matrix of the dynamical modes; 
the remaining terms, denoted by the dots, are omitted for brevity. Specifically, the terms included in (|69p are those 
proportional to the coefficients a.y , 5^ , and Cy in eq. (jl9p , where the quadratic action is written at a formal level. These 
arc the only terms entering in the kinetic matrix K of the dynamical modes, see eqs. (1231) and (|24p . It is immediate 
to compute K from ((69|) . However, the explicit entries of this matrix are very involved, and not illuminating. For 
this reason, we do not report them here. We can however compute the eigenvalues of this matrix numerically for any 
given choice of parameters (namely, for any background evolution, and momentum of the perturbation) . In addition, 
it is possible to obtain a simple approximation for the determinant of K 

"^^^ah^) ~ 192 p^ 32p2p6 ^^\pi) ^"^> 

where Hq has been defined in eq. (|68p . This expression is accurate in the sub-horizon regime during inflation, in the 
limit of small anisotropy. More specifically, we obtained it using the same steps outlined after the analytic expression 
for the determinant of the model considered in the previous Subsection, eq. ([64]) . It is worth noting that (jTO]) differs 
from ([64|l only by an extra factor of 1/2; this is the original factor in the kinetic term of the infiaton (which is the 
only additional field in the model we are considering in this Subsection). This suggests that the perturbation Scj) is 
decoupled from the other ones at leading order. 

In figure [5] we present the results of a numerical evolution for a given set of parameters and initial conditions 
(starting from inflation, in the slow roll regime). The left panel shows the evolution of various background quantities. 
The right panel shows instead the determinant of the kinetic matrix for a mode with = lOOiJ, px = 200H at the 
initial time (both the exact value, obtained from a numerical evaluation of the kinetic matrix, and the approximated 
value (|70p . are shown). We see that the determinant vanishes and becomes negative in the sub- horizon regime, when 
H/p - B . 

In Figure [H] we show the evolution of the four eigenvalues of the kinetic matrix; notice that there is an additional 
dynamical mode, supported by the inflaton field, with respect to the model studied in the previous Subsection. 



23 




5 10 15 20 25 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 

mt H/p 

FIG. 5: Results from a numerical simulation starting at f = from <^ = 16 (providing about 60 e- folds of inflationary expansion), 
Bi — 0.1 (providing a ~ ~ 1% anisotropy). More precisely, we have considered a massive inflaton potential, with the inflaton 
mass equal to m . Left panel: inflaton (in units of Mp), hubble parameters (in units of m), and dimensionless rescaled vectorBi. 
The anisotropic rate h and the vector are rescaled so that they are visible in the figure. Right panel: determinant of the kinetic 
matrix of the perturbations, for modes with initial momenta PL,in = lOOi/in, PT,in = 200-H'in, ■ The red curve shows the exact 
determinant, while the green points show the approximate expression given in eq. ((70}. The determinant vanishes at the time 
mt ~ 0.16. 



However, as for that case, one mode is a ghost for some time. We know that the system of hnearized equations for 
the perturbations become singular when the eigenvalue Ai crosses zero. We expect that also the solutions diverge at 
that point, analogously to the study of the previous Subsection. 



D. Vector inflation 



We now study the simplest realization of the idea of vector inflation proposed in It is characterized by three 
mutually orthogonal vector fields nonminimally coupled to the curvature. The three fields have equal vev, providing 
a FRW background. Despite the background solution is isotropic, the model suffers of the same ghost instability as 
the models studied above. The discussion is divided in two parts. We first introduce the model, and discuss the 
background evolution. We then study the spectrum of perturbations around this background solution. 



1. The model and the background solution 
The action of the model is given in cq. (fT^ , with ^ = — 1/6 , and = 3 : 



S = I d^xV—gfl 



_1 p(a) p{a)fii' _ 1 
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(71) 



The background vev of the three vectors is chosen as in eq. (jl3p . while the background geometry is ds^ 
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FIG. 6: Eigenvalues of the kinetic matrix, for the same choice of parameters and initial conditions as in the right panel of 
Figure [5] The eigenvalues A2,3,4 are always positive, so that the sign and the behavior of the determinant are determined by 
Ai. This eigenvalue is initially positive, crosses zero at H/p ~ 0.082 (see the previous figure), and diverges at H/p ~ 4.6. In 
the y axis, we have used a linear scale inside the interval [—0.01, 0.01], and a logarithmic scale outside. 



-dt^ + (t) dx^ . The system is governed by the equations 
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(72) 



For the background under considerations, the 00 Einstein equation and the i — th spatial component of the equation 
for the i — th vector give 



B + 3HB + B = 



(73) 



As always for a FRW geometry, the diagonal ii Einstein equations does not provide additional information (due to a 
nontrivial Bianchi identity). The remaining equations vanish identically. Upon the identification B = (pctf/ (VSAfp), 
we recover the same equations as those of chaotic inflation driven by a minimally coupled scalar field 0off (where 
the suffix stands for "effective") with potential V = m^(/)^jj/2. Inflation is therefore characterized by the slow roll 
evolution 
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(74) 
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which apphes as long as the slow roll conditions 



< 1 



(75) 

are valid. After that, B performs damped coherent oscillations around B = (see for instance p^): 

S= 1:5^ + sin (mt + eo) (76) 



3mt \ m2 t 

where is a phase (irrelevant for this discussion). The coherent oscillations provide a stage of effective matter 
domination {w — once averaged over a few oscillations) [i^. It is expected that the vectors then decay into the 
visible matter (either perturbatively, or nonperturbatively) giving rise to the radiation dominated stage of cosmology. 

From the slow roll conditions, we see that ^ during most of the inflationary stage. At these times, the 
mass term for the vectors is negative, —R/6 + m? ~ —2H^ + < . However, R decreases as B rolls towards the 
origin, while m remains constant. We find numerically that —R/6 + ni^ = when B ~ 1.048 . This happens towards 
the end of inflation. 



2. Ghost instability from the quadratic action 

We now study the perturbations of the model l(7T|) around the background solution just discussed. Since the 
background geometry is FRW, one may choose to adopt the standard decomposition of metric perturbations, and 
decompose them into scalar, vector, and tensor modes with respect to 3d spatial rotations [s?!. While this is possible, 
one should however bear in mind that, contrary to the standard case, these three groups of modes cannot be studied 
separately even at the linearized level. Consider for instance the tensor perturbations hj,^ , introduced as the traceless 



(/i^'^ = 0), and transverse {dj hf-^ = 0) part of the spatial metric perturbations, 6gij = hf^ [37[. While in the 



case of scalar field inflation these modes are decoupled from the other perturbations, and can be studied separately, in 
the case under consideration they are coupled to the perturbations of the vector fields. Namely, we find the following 
coupling in the quadratic action S2S of the perturbations: 



J d^x a^[(^B + HB) 6Af + (2H^ + H- m^) B 5A^ + {t ^ j)] hj^ (77) 



As a consequence, the perturbations of the vector fields must be included in the linearized equations of motion for 
the tensor metric perturbations. In turns, the perturbations of the vector fields are coupled to the scalar and vector 
perturbations of the metric. One then finds that all the perturbations of the system need to be studied together, even 
at the linearized level. 2^ 

This makes the study extremely hard, and indeed no complete computation exists up to date. The metric has 
10 perturbations, while each vector introduces 4 perturbations. Thus, in this simplest realization of vector inflation, 
one starts with a system of 22 coupled modes. Four perturbations can be removed by fixing the freedom of general 
coordinate transformations, leading to a system of 18 coupled modes. Not all these modes are dynamical. The 
nondynamical modes originate from the initial perturbations Sgofi and SAI^\ which enter in the quadratic action of 
the perturbations without time derivatives. As we now show, it is possible to choose gauge invariant perturbations 
which associate a gauge invariant combination to each of the Sgofi and SA'^^ modes. These seven gauge invariant 
combinations are also nondynamical fields. We integrate them out as outlined in Section [TTTl and obtain the quadratic 



See I41II for more detailed studies of initial and slow-roll conditions for vector inflation. 

Ref. [32l | studied the tensor modes alone, arguing that the effect of the coupling to the ^A^"^ perturbations can be "averaged away" in 
the limit of many vectors. However, each vector introduces a coupling of hTT with its own perturbations, and perturbations of different 
vectors cannot cancel against each other in the study of the spectrum (and, therefore, of the stability) of the theory. 
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action for the dynamical modes. From the study of the kinetic matrix of this system, we find that three modes become 
ghosts for some time during inflation. 

Since the usual decomposition in scalar/ vector/tensor modes does not provide decoupled sets of perturbations, we 
do not employ it, and simply decompose the perturbations as 
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where A^^^"'^ denotes background values of the vector fields. The paramctrization of the ^1733 component has chosen 
so that 6gl — —2 T . Moreover, for algebraic convenience, some Sg^i, entry has been written with spatial derivatives 
(in practice, we assume that the momentum of the modes - after Fourier transforming - has nonvanishing components 
in all three directions, k^, ky, k^ 7^ ). 

We need to fix the gauge freedom associated with general coordinate transformations. Consider the infinitesimal 
transformation x*^ + , under which the perturbations of the metric and the vector field transform as Sg^^i, — > 

'^5^11^ ~ g'iS},a^°' — ffla'?^ — J ^'^d (5A|f ^ — + ^ — — ^ ■ Due to the assumption of the modes we 

are studying, we need to consider infinitesimal transformations with nontrivial spatial dependence along all the three 
directions. We can therefore parametrize the transformation parameter as — (^*', 82^^, 83^,^) . The explicit 

transformations of the modes in (l78l) are then 
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(81) 



are gauge invariant (all these modes have mass dimension +1 ). 

We expanded the metric and vector fields according to ((78|) . and computed the quadratic action for the perturbations. 
We verified that the perturbations rearrange so that the quadratic action can be written solely in terms of the gauge 
invariant modes (|8ip (this is a nontrivial check on our algebra). In this way, we have eliminated the redundancy 
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associated to general coordinate transformations. The action in Fourier space reads 
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where we have used the physical momenta pi = ki/a. Eq. (j82p actually is not the full quadratic action of the 
perturbations, but some terms (denoted by dots) are omitted. Let us clarify this. We find that no time derivatives 
of the combinations Xi, dp"^ enter in the quadratic action (neither in the terms shown here, nor in those omitted). 
These are the nondynamical gauge invariant modes of the system. The remaining modes are dynamical. Eq. (j82p 
is the action for all the gauge invariant modes of the system (both the dynamical, and the nondynamical ones). 
However, it contains only the terms that contribute to the kinetic matrix of the action for the dynamical modes, once 
the nondynamical modes have been integrated out. These terms are given without any omission, nor approximation, 
so that eq. (|82p contains all the necessary information for the exact computation of the kinetic matrix of the dynamical 
modes. 

Before integrating out the nondynamical modes, we can easily see that 5 dynamical modes decouple from the 
remaining ones in the part of the action shown. There are 9 dynamical modes in the spatial perturbations of the 
three vector fields, d^"'' . These perturbations enter in the action ([5^ with a diagonal quadratic term (the second 
term, cx Id'"''!^). Then, they are coupled with the remaining modes only in the third line of ([5^ ). We see that, out 



of the nine modes d^"'' , only the four linear combinations Pjd^"'' (a = 1, 2, 3) and 0.^' are involved in these couplings. 
The remaining 5 linear combinations are decoupled. Therefore, we can rotate the fields d^"'' into the coupled and 



(a) 



Notice that the procedure we adopted is equivalent to choose the gauge E\ = E2 = E3 = T = 0. It is easy to see from I I79I I that this 
choice can always be made, and it completely fixes the gauge freedom. 

Notice that they corresponds to the nondynamical perturbations Sgof^ and iJAq"' in the gauge E\ = E2 = -B3 = T = 0. 
For clarity, compare with the formal discussion of Section IlIII The action I I82I I given here corresponds to the formal action I I19I I. with 
only the terms proportional to the coefficients aij, bij, and Cij included. Those are the only terms necessary to compute the kinetic 
matrix K of the dynamical modes, cf. eqs. II23II and II24II . 
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decoupled linear combinations: 
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where for brevity we have defined pij = ^pf + p^ 

The combinations vi, . . . , £4 are the coupled modes, while V5, . . . , vg are decoupled from the remaining perturba- 
tions. One can check that the matrix relating the vector la^^^"* , dj^'' , dg^^ , a^i \ . . . , a'^'^ | to the vector {vi , . . . , wg} is 

orthogonal. Therefore, the quadratic term oc ^ |q;^°''P becomes oc \vi\'^- As a consequence, the modes v^, . . . , vg 
enter only through the term 

and are therefore decoupled from the other modes in the part of the action shown in ([5^ (we stress that this part is 
all we need to compute the kinetic term for the dynamical modes). Since the coefficient in front of each term in the 
sum is positive, we see that these modes are well behaved (positive energy) excitations. Let us now focus on the 
coupled system of the remaining modes (which includes the remaining 6 dynamical modes, and the 7 nondynamical 
modes). The next step is to integrate out the nondynamical modes as outlined formally in equations form (|19p to (j24p . 
Doing so, we end up with a rather long kinetic matrix for the six dynamical modes j^*, S, iii} {i — 1, 2, 3, 4). We were 
unable to perform analytical simplifications of this matrix, as those leading to the expressions (I64p and ([70]) for the 
simpler models studied above. However, we can evaluate this matrix numerically for any given choice of parameters. 

More specifically, we show the results of a numerical evolution starting from Bin = 10, i?in = according to 

the slow roll conditions ([Ti]) . and with momentum of the modes initially satisfying pi — 100 H, p2 = 80 H, p^ — QOH 
(these values have be chosen only for illustrative purposes and have no particular relevance, other that we require 
the mode to be well inside the horizon initially; we verified that other choices of the momenta lead to the same 
conclusions as those shown here) . We show in Figures [7] and [5] the evolution of the eigenvalues of the kinetic matrix 
(the second Figure is a close-up of the first one where the eigenvalue A3 crosses zero) close to horizon crossing. We 
see that two eigenvalues Ai_2 are initially negative; the corresponding eigenvalues are ghosts at these times. These 
eigenvalues become positive without crossing zero (they diverge). A third eigenvalue A3 is initially positive, but it 
becomes negative crossing zero. It then diverges and becomes again positive. As the formal eqs. ((22|) show, the 
linearized system of equations becomes singular when this eigenvalue crosses zero. We expect that also the linearized 
solutions diverge at this point (cf. the explicit solutions given in the model of Subsection IV B[) . 



The overall factor appearing in l|82| l has not been included in this computation. This is irrelevant for our discussion, since this factor 
simply rescales all eigenvalues by a common positive number. 
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FIG. 7: Eigenvalues of the kinetic matrix for the dynamical and gauge invariant perturbations 5', E, vi,...a for the model 
of vector inflation (I7ip . and for one specific choice of initial conditions (see the main text). Since H/p increases with time 
during the stage shown, we use it as a "time variable" in this and the next Figure. The eigenmodes corresponding to Ai,2 are 
ghosts at the earliest times shown (low H/p). The mode corresponding to A3 becomes a ghost close to horizon crossing. This 
eigenvalues (and the determinant of the kinetic matrix) crosses zero at this point, signaling an instability of the system also at 
the linearized level. The kinetic matrix, and its eigenvalues, are dimensionless. Notice that we use linear units in the interval 
[—0.0001, 0.0001], and logarithmic units outside. 

The study so far concentrated on the nature of the modes at horizon crossing, in an inflationary regime for which 
the total mass term of the vectors was negative, — i?/6 + ~ —2H^ + m? < 0. However, as we discussed after 
eq. ([75]) . —R/6 + m? vanishes towards the end of inflation. We studied the behavior of the eigenvalues of the kinetic 
matrix also around this point. We considered the same background evolution as for the previous plot, but smaller 
values of the momenta, so that H/p is not exponentially small at the times shown (we want to avoid that our results 
are affected by numerical inaccuracies). As shown in Figure [9l we find that two eigenvalues cross zero precisely when 
the total mass vanishes. Also at this point, the system of equations for the eigenvalues becomes singular We expect 
that the linearized solutions diverge also at this point. 

It is interesting to compare the behavior of the eigenvalues shown in these Figures with that obtained for the 
previous models. We find that the two eigenvalues A1.2 behave precisely as in the case of zero vev studied in the 
previous Section: they are negative in the deep sub-horizon regime, they diverge close to horizon crossing, and they 
cross zero when the total mass vanishes. On the contrary, the eigenvalue A3 behaves precisely as in the cases of a 
single vector with nonvanishing vev studied above (cf. Figures [3] and [6]) : it is positive in the deep sub-horizon regime, 
it crosses zero close to horizon crossing, and it remains negative for some time afterwards. It appears from these 
behaviors that the mixing with gravity affects only one linear combination of the three ghosts. 

We conclude the present discussion with some remarks on the previous study |33| of perturbations of vector inflation. 
We already discussed in the Introduction while the computations of ^33| - being limited to the linearized equations 
for the perturbations in either the sub-horizon or the super-horizon regime - cannot show the ghost instabilities found 
here. Here we want to reply to some specific comments on our previous works contained in [3^. Most of these 

remarks are answered by the fact that, in our previous works, we only provided arguments for the presence of ghosts in 
the model of vector inflation, deferring the explicit computation to the present work. It is mentioned in [s^ that the 
ghost may be an artifact of having expanded the vector field in transverse and longitudinal part according to Aj^ + dfj^cj). 
This would introduce additional time derivatives, which could affect our findings. This Stuckelberg decomposition 
was introduced in [2^, |3l[ only as the simplest way to elucidate the problem. The computations presented here do not 
introduce additional time derivatives in the parametrization of the perturbations, and lead to the same conclusions as 
the much simpler Stuckelberg analysis. It was also pointed out in ^] that the complete computation, with gravity 
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FIG. 8: Close up of the previous Figure where A3 vanishes. A further close up shows that A5 and Ae, and, later, A3 and A5 do 
not cross each other, although they appear to do so in the Figure. 



perturbations also included, was in order. This is precisely what is performed here. Another comment of [33| was 
that the model of vector inflation should only be regarded as an effective field theory, valid only at small energies. 
The (unknown) UV completion of this theory may be without ghosts. We agree with this claim. ^® This was precisely 
the point raised in our previous works. As we have seen here, the ghost instabilities appear during most of the sub- 
horizon regime, and close to horizon crossing. Any UV completion needs to be relevant at these stages. Therefore, 
the effective theory of vector inflation cannot be used for the study of cosmological perturbations in the sub-horizon 
regime. Since this regime is crucial to obtain phenomenological results, this invalidates any prediction of the effective 
model [2^. Finally, it was argued in [ssj that the instability may simply be due to the growth of the scale factors, 
and to a wrong rescaling of fields. This is not the case, since the linearized system of equations - and, most likely, its 
solutions - diverges at some finite moments of time, while the scale factor remains finite. 



VI. DISCUSSION 



Although the paradigm of inflation is well supported by the observational data, we still do not know the actual 
particle physics mechanism behind the inflationary expansion, and the generation of cosmological perturbations, ft 
is customary to parametrize our ignorance in terms of scalar fields: they may be fundamental fields, or simply order 
parameters which provide an effective description of some degrees of freedom in the theory (for instance, the brane- 
antibrane separation in some string motivated models of inflation). However, it may well be possible that these two key 
elements of cosmology are due to higher spin fields. Vector fields are probably the simplest possibility after scalars. 
They can in principle leave a distinct signature from the scalar case, since a nonvanishing spatial vev of a vector 
breaks the isotropy of space. This can provide anisotropic expansion, and / or generate a spectrum of primordial 
perturbations that breaks statistical isotropy. The main obstacle faced by explicit realizations of this idea is that, 
in the standard case, vector fields decrease too quickly due to the expansion of the universe. Therefore, suitable 
modifications of the standard action need to be made. 

Recently, a class of models was considered in which the vector is nonminimally coupled to the curvature, C D 



Apart from the fact that we are no aware of any well behaved UV completion of a theory with ghosts. 
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FIG. 9: Two eigenvalues vanish when the mass term = — R/Q of the vector fields vanishes. The mass is shown 
in units of m^/300 . We have B ~ 1.048 when the total mass vanishes (this occurs still during infiation). The mode has been 
chosen to be outside the horizon when the total mass vanishes {H/p ~ 1872 at the time shown, and fci : fc2 : fcs = 10 : 8 : 4). 

R/12 A^^ . Indeed, for this specific coupling, the vev of the vector evolves as that of a minimally coupled scalar 
field; this offers the possibility of realizing an inflationary background, with a controllable anisotropy [1^, [1^. In 
addition, the transverse perturbations of the vector behave as the perturbations of a minimally coupled scalar; this is 
the basis of the mechanism of vector curvaton (30,] for the generation of a nearly scale invariant spectrum of primordial 
perturbations. To our knowledge, the analogy between the R/\2 A^ coupling and the minimally couple scalar field 
first appeared in the 1987 work by Turner and Widrow [2?], where it was exploited for the generation of a primordial 
magnetic field during inflation. The renewed interest in this mechanism is mostly motivated by some features in the 
WMAP data that hint for a small break of statistical isotropy. 

All of the mentioned works suggested new interesting mechanisms, and a complete check of the stability of these 
proposals was beyond their scope. It is tempting to assume that these models should be well behaved due to the 
strong analogy with the minimally coupled scalar fleld case. There is however a crucial difference between these two 
cases, and between the case of a minimally vs a nonminimally couple vector; it is due to the longitudinal vector 
polarization. In the above works, the vector has a U(l) invariant kinetic term C D —l/AF^uF^'^ . If only this term 
was present, the vector would only have the two transverse polarizations. However, the nonminimal coupling to the 
curvature breaks the U(l) symmetry, and gives rise to an additional, longitudinal, polarization. The nature of this 
mode is controlled by the sign of the mass term, which, for these mechanisms to work, needs to be negative. For the 
scalar case, a negative mass squared means that the fleld is a tachyon; for a vector fleld, a negative mass squared 
means that the longitudinal polarization is a ghost. 

Motivated by this initial consideration, we studied the stability of this class of theories. We did so by computing 
the free action for the dynamical (physically propagating) modes of such theories, around the background solutions 
considered in the various proposals. The sign of the eigenvalues of the kinetic matrix of these action indicates whether 
the corresponding eigenmode is a positive or negative energy excitation. Our computations conflrmed that there 
is a ghost for each nonminimally coupled vector fleld in the model. As we already mentioned in the Introduction, 
theories with ghosts are consistent only as effective theories, valid below some energy scale A. Inflationary predictions 
heavily rely on the initial conditions; for instance, the vector curvaton mechanism of [30| results in a scale invariant 
spectrum because of the speciflc choice of initial adiabatic vacuum. This choice is made in the quantized theory for 
the perturbations, which is performed in the deep UV regime (energy ^ H~^) [s^l- In presence of a cut-off, the initial 
adiabatic vacuum cannot be imposed at arbitrarily early times, and, depending on the precise numerical value of the 
cut-off, it may not be possible to impose it at all. This casts doubts on the phenomenological predictions obtained 
for such models. 
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In fact, all theories with an explicit mass M for the vector require a cut-off which makes them invalid at high 
energies, irrespectively of the sign of the mass term. We can see this based on the behavior of massive vector fields 
at high energies. The explicit mass breaks the gauge invariance in a hard way. It is well known that, in such case, 
the interactions of the longitudinal bosons violate unitarity at a scale which is parametrically set by M, leading to a 
quantum theory out of control. For the present models, M is the Hubble rate or below, so that the entire sub-horizon 
regime may be ill-defined. At high energies, the longitudinal mode will also interact strongly with the other fields 
in the theory (this will renormalize the coupling constants). Then, depending on exactly when this happens, the 
quantum theory of the perturbations may be out of control throughout the entire short wavelength regime. If this 
is the case, all initial conditions would become unjustified, and the theory would lose its predictive power. Although 
this problem is present for both positive and negative mass terms, a theory which has a hard vector mass and a ghost 
is more problematic than a theory with only a hard vector mass. The most immediate UV completion of a theory 
with a hard mass is through a higgs mechanism. The mass would be then due to the vev of a scalar field that becomes 
dynamical above the scale M . In this way the theory remains under control also in the short wavelength regime, and 
one can apply all the standard computations valid for scalar fields during inflation. However, if needs to have the 
wrong sign, the scalar field in this UV completed theory needs to be a ghost. In fact, we are not aware of any well 
behaved UV completion of a theory with a ghost. 

We stress that this instability is not related to the classical behavior of the solutions of the linearized equations for 
the perturbations, and it would be present even if the latter remained finite. However, we argued that, for the models 
considered here, also the linearized solutions diverge. This is a second instability that adds up to the one we have just 
discussed. This instability appears because the eigenvalues corresponding to the ghosts do not remain negative over 
the whole evolution, but change sign, and cross zero at some finite moment of time (there are two such moments 
in the model of vector inflation considered here). We showed that the linearized equations are singular at t*, and we 
expect that the linearized solutions diverge for t — > . While it is possible that the divergency does not occur at the 
full nonlinear level, this instability also invalidates all the phenomenological signatures of these models which are based 
on the linearized computation (as for instance the primordial spectrum of perturbations). We solved the linearized 
equations only in the simplest cases of a vector field with no vev, and of a vector field with vev plus a cosmological 
constant. We did not solve them for the models i29j|. We have shown however that the linearized equations become 
singular also in these cases. It is important to stress that even if, due to some unexpected cancellation, the solutions 
to these equations would remain finite, this would not eliminate the ghost instability that we have discussed in the 
two previous paragraphs, and that we have proven to be present for all the models studied here. 

We conclude the Discussion with some remarks on models different from those considered here, and in our previous 
works [25l[3]|. The instability we pointed out motivates the study of such models, as for instance vector fields with 
nonstandard kinetic terms [ly, |4^|4^ (although, some of these proposals are also unstable), nonabelian vectors with 
nonminimal coupling to the curvature [4^ spinors [i^l, or higher p— forms [2l|,[4l, [4^ Of particular interest, 
in our opinion, is a class of models characterized by a function of a scalar field multiplying the kinetic terms of the 
vectors, I (</>) F^^" F^i, , but no potential term for the vector [i^,!!^]- U(l) invariance is preserved in these models, and 
the problematic longitudinal mode is absent. A complete study of the cosmological perturbations (conducted along 
the lines described here) is the next necessary step for obtaining the phenomenological predictions of these models. 
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The stability analysis performed here applies also to the nonabelian model of [46l | if the vectors have no vev (since, in this case, the 
nonabelian structure does not affect the quadratic action for the perturbations); however, the case with a nonvanishing vev requires a 
separate investigation. 

The vector case is p = 1; Ref. figjl generalized the arguments of [25ll3l|. and pointed out that also the p = 2 case contains ghosts. 
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APPENDIX A: DETAILS OF THE LINEARIZED COMPUTATION OF SUBSECTION IVBl 



This appendix contains the details of the stabihty analysis of the model of Subsection IVBl The explicit forms of 
the linearized equations (|53p are 
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where the index i on the equations spans over the i = 2,3 isotropic spatial directions, and where the momenta are 
defined in eqs. ([5^ and ((55|) . 

Only the 2d scalar perturbations are included in this computation. More in general, if we include both the 2d 
scalar and 2d vector modes, the perturbed equations carrying i or j indices split in two separate parts, one for the 2d 
scalar, and one for the 2d vector modes. Although we have not written the 2d vector parts of the above equations, we 
have explicitly checked that they are decoupled from the 2d scalar parts. The 2d vector parts are not related to the 
instability we have demonstrated in the main text; therefore we do not discuss them here. Another property to be 
noted in the system (jAl[) is that not all the equations are independent. Using the perturbed Bianchi identities, it is 
possible to obtain eqs. Eq^, Eqj^j from the remaining ones in (|A1[) . This remaining equations are the set of linearized 
equations ([55]) that we have chosen to solve in the main text. 

For brevity, we have grouped some of the extended terms that depend on the background quantities, which appear 
both in the action and in the linearized equations. We have denoted them with calligraphic letters D, Ai. These 



•^^ One could have chosen a different but equivalent subset of independent equations. Our choice is related to the fact that the linearized 
equations of motion obtained by extremizing the quadratic action II63I I precisely gives the set of equations I I56I I. 
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1 



18 + 3Bf + 2BI 



-1 



18 + Wl + 2Bf 



3p. 



l^^ + }^h'-^hH] (3m'-6H'-^h^ + 15hH) Bf-^hB.B^ 

2 8 2/ V 4 J 2 

+ ^{6h- 7H) Bi Bf + 27 (8/i - 3H) Bi Bi + 27 (^3H^ - 3h'^ - ^ Bf 
+3 (em^+'-^h^-'-lH^ + fBf^ Bf 

f - 9. h) Bl ^ 27 (2. - ^) + 3 - 5.) - ^ (6. + H) B. Bl 

+9 Um'-&h' + l2hH-^H^+'^Bi\ B, + 3 f^m^ + ^h^-UhH-^H^+^-Bl] Bf 



15 



pI (18 + 3Bf + 2Bf) 



Pt 

-1 

pI (18 + 3Bf + 2Bf) 



— /i^ - 3/iff Bj' - UhBiBi + {2H-7h) B^B^ 



+9 (m^-2h^ + 4hH-3H^ + -Bf] + (-m^ + —h^-12hH--H^ + - Bf) Bf 
\ 2 J \2 2 2 4 / 

{^h'^ -2lhH + 3H^^ Bf-36hBiBf + 3 {2H -7h) BiBf 
+9 (3m' + 6h' -QH' + l bA Bf + 3 (^m'' + ^ h'' - 18hH - 3H' + ^ Bf\ Bf 



(A2) 



18 + 3^2 + 2Bf 



[m 



+ 5/i^ - AhH) Bf -AhBfBi-2 {AH - 19/i) Bf Bi +24{H- 2h) Bi Bi + 3Bf 



+ I 



91, 



2 ... ^h'^+46hH- 12H'' + Bfj B{ - [3m^ + 2lK^ + 12hH- 36H^ + ^Bf ] Bi 
{2m' + - AhH) Bf -3Q{H- 2h) Bi + 4{H - h) Bf Bi + 6 {3H - 16h) Bf B2 



18 + 3Bf + 2Bf 



2pi (18 + 3Bf + 2Bf) 



-9 (m' -13h' + 12hH -2H' -^Bfj Bf + 18 (m' + h' + 4h H - 6H' + '^Bf ) Bi 



- 3 (4iJ - 5/1) hBl+A{H- 5h) Bf Bi -36{H - 2h) Bi - 6 (6/1 + H) Bf Bi 



+36 (m' -2h' + AhH - 3H' + ^Bf) Bi+Q fm' + 13h' -8hH-3H'+ '^-Bf ) Bf 



36 



Mao, 



M 



18 + H 


h2B4 


2 




18 + iBl H 


h2B4 


12 




18 + iBl H 


h2B4 


4 





lO/i^ - Sh H) Bf -ShBfBi-e (2H - 5h) Bi Bi + 9 [ - 5h^ + Ah H - -Bf 



-3 ('^m^ 



13, 



IQh H - 6H^ + Bl] B't 



(m^ + 22/1^ - Uh H) Bf - 8/1 Bi - 6 {2H - 5/i) Bi Bi + 9 m 



7h^ ~2hH- -B 
3 



+3 



7hH -GH"^ + B^] Bf 



Ah H) Bf -AhBlBi + 2 {7h - AH) Bf Bi + 6 {7h - 2H) Bi Bi - WBf 



15m^ 



5 • . 



11 



57h^ + 60hH- mn'^ -7;Bi] B'f + [ - ^h^ + 10 hH - 12H^ + B^ 



18 + 3Bf + 2B1 



2 TO^ 



7, 



2hH] Bl + -{h~ 6H) BfBi+9 {7h - 2H) Bi + {2H - llh) Bf Bi 



1 



+36 \^m^ - -h^ + AhH- -H^ + -Bf j Bi + 9 |^m^ + 2/i^ - 2i/^ + 9 ^1 j ^1 

As we have discussed in the main text, we proceed by solving the second of (j56p for the mode x, and inserting the 
solution back into the rest of the equations. Next, we differentiate equations EqQQ, EqQj, Eqp (with the solution for x 
given in (|57p inserted in them) in order to obtain first order differential equations for the modes $, B, dg. Combined 
with Eqj^^i Eqj^, Eq^ (again with the solution ([57]) inserted in them), these equations form the set of equations to be 
solved numerically, summarized in matrix form in (j59p . Here we give the detailed expressions of the terms appearing 
in this system of equations. 

We first define the following useful combinations of background dependent terms which frequently appear in the 
linearized system: 



V = {6 + Bf) pI^^pIV^ 
Bi Bi 



n 



H + h 



Q + Bl 



(A4) 



where P^^^, is defined in (jA2p . The time derivative of T) and Ti. are also useful, which is explicitly given by 



V ^ -2p| (6 + Bl) {H + h) + Wpl {H - 2h) V^^ + 2p| Bi Bi - 8pl U 



XX 



n 



BiBi 
6 + Bf 



n + H + h + 2{H + h) 



Bl Bl Bf + Bl Bl 
e + Bf ^ 6 + Bf 



(A5) 
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where Bi, H, h are obtained from ([49]): 



H = 



h = 



Bi = 



2 


18 H 




+ 2BI 




1 




18 H 


-3^2 


+ 2Bf 




1 




18 H 




+ 2Bf 



-^h^ + 7hH -3hA B^ + 12hBiBi + {7h-2H) B'f Bi - ^ (eh^ + Bf 



(l8h^-12hH+'^Bl)Bf 



3 (2iJ - 5/i) hBt + 8hBfBi+6 {2H - 5/i) Bi Bi + 6 [Bf - 9h H 



11 



-6 ( -2h^ + —hH - 3H^ + -Bf ] Bf 



-2{m^ + lOh^ -8hH) Bf + 2 {8h - 3H) Bf Bi - bAH Bi + 15{H - 4/i) Bf Bi 



-18 ( - + 4:hH~ -Bf ] Bi - 15 [ m 



5 5 5 



We also need the explicit expressions for t^xx ^^'^ T^aoao- 

3 + ABf 



H -2h- 



18 + iBf + 2BI 



4 Bi Bi 



XX 



9 



Vi im + 3Bf + 2Br, 



2 ( - 7hH + H^j BlBi-12hBfBl + ^{2H-7h) Bf Bf 



+6{m-' + 2h^ - 2H^ + -Bf \ Bi Bi + 4 , ^ 



^Bf] BfB, 



2 ■ 7 : 



31. ; 7 



AhBfBi+[-H--h] BlBi+\—hh-- [Hh + hH] + -HH ] Bf-AhBfBi 



6 



7. 



-H - -h] B° Bi + [12 [hh- H H] +3BiBi) Bf 



+ [17hh-6(Hh + hH] -2H H +^BiBi] Bf 



i'aoaa — y^XX 2^^ ^^XX 



(A7) 



Now we give the exphcit form of the matrix Ai^ used in ([5; 



2 , 2piBf 



Kl2 
Kl4 

Kl6 



I? 

Bi 2pi Bl 
2 V 

Bf\ , 1 ■ 2pj Bf (6 + Bf) 
3+^j h+{3-Bf) H + -BiBi- ' ^ '-^n 

3 Bf plBf {6 + Bf) 
2 4 

V 



2V 



38 



H21 
K25 



Bi ^ ^pjBf 
3 



31? 



= 1 



3V 



K26 = 1 



K41 = 1 + 



6 



h 



1 



Bf 



3V 



H+-B,Bi 
6 



2pi 5^(6 + 5?; 



31? 



n 



2 31? 



K44 — 3-ff 



1 ^. 



K45 



Bl 
6 



j /i^ - - Bf + 2 Bl [Bl H + Bl j h 

. I ■ p] {Q + BlY 
(2g-/z) + -BiBi + ^^ 32? ^ 



2pi (6 + Bf ) 
31? 



K46 = Bl - — r- 



H\ Bl 4p2Bi(6 + Bf; 



H2?a 



K51 



K54 



1 - b; 



31? 

plBl {6 + Bf 



1 / ^ pI (6 + Bf 



6 



I? 



61? 

Bl 



1^56 — 



^ (6 + Bf) (2B-;i) + iBiBi 

pi{(^ + B!) ^ 

p| Bl (6 + Bf ^ 



pI {(^ + Bfy 

6V 



n 



3V 



I?n 



K61 



K63 



'^65 



31? 

L 

32? 



I?> 



K62 — 1 + TTTT ^XX 



K64 = Bl (B - 2/i) + Bl 
2pi (6 + Bf^ 



3BiI? 



I? 



8pi (6 + ^?) ^^ 
3BiP 



XX 



^ ^ p| 2p^(6 + Bf; 



Pi 



31? 



I?n 



(AS) 



Finally, we explicitly write down the right hand side of ([59|l involving the functions fi where i = 1 . . . 6. Each fi is a 
linear combination of the variables 



Fi = di, d, di, d, l>, B, do} 



(A9) 
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The explicit forms of the functions /i, . . . /e can then be expressed as 



(AlO) 



The coefficients Aij depend entirely on the background quantities. They are given by 



V 
V 



H 



Ah Bi 



V 



Bi V 



^13 = 

Au = Pt^i - + 



V 



I, 



V 



H -2h-2 



Bi 
'Bi 
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Bi 
Bl B, 
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Bl 



Bl Bl 



B-2h-^]+iH-2h)^^-^^ 
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17 



r I _ + 3/j2 _ ^ j ^ 1 (9^2 _ 2pl-p^^ + 15/i2 - 12hH) Bf - 3hBi Bi 



M2 
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2plBf {6 + Bl) 



V 
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A 
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A28 
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4b? / Bl B? 



\ 
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Bl V j V 
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A2 



3D 
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H Bi\ ■ H , Bl Bf Bl 
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A26 = Pt 
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3 ; 32? 
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H + n 
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3D 



2V 2Bi 2Vo 
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^31 = ^32 = A34 = A37 = 

A33 = -3{H- 2h) 

A35 = pi 

A3Q = -pi - Maa 

2 
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2 ^. 
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(All) 



Thus, we have the fuU form of the equation (j59p . which we integrate numericaUy. 



APPENDIX B: EARLY TIME CANONICAL ACTION AND INITIAL CONDITIONS 

We discuss here how we set the initial conditions for the perturbations entering in the hnearized system (j59p . As 
in the standard case [svt , the initial conditions follow from the quantization of the quadratic action for the dynamical 
modes. As we show in Section Hill the quadratic action - formally written in eq. ([M]) - is obtained by integrating the 
nondynamical fields out of the quadratic action for the perturbations - formally written in eq. pop . 

The quadratic action of the 2d scalar perturbations of the model (|46p is given in eq. (|63p . The corresponding action 
for the dynamical modes is extremely lengthy, and we do not explicitly write it here. Fortunately, to set the initial 
conditions for the perturbations we only need the leading terms in an early time expansion of this action. We are 
interested in the phenomenologically relevant case of moderate anisotropy {Bi <C 1), for which Ha — Hb are nearly 
constant. Therefore, the two scale factors a and b grow nearly exponentially with time, and p/H (where p is either 
the longitudinal or the transverse component, and H is either Ha or Hb) is exponentially large in the asymptotic past. 
Therefore, the early time expansion of the action coincides with the sub-horizon p/H 3> 1 expansion, exactly as in 
the standard case. 

Specifically, we first compute the exact matrices K, A, fl^ that form the action for the dynamical modes (cf. eq. 
(1221)) J and then expand them for p/H 3> 1. Since the resulting expressions are still quite involved, we further expand 
them for i?i <C 1. This procedure is legitimate provided that 

— < Si < 1 (Bl) 
P 

which, as we remarked, is always true in the case of small anisotropy, and for sufficiently early times. At leading 
order, we obtain the expression 

5ca„^^ J d^kdt [\H+\^-p^\H+\^ + \A+\^-p^\A+\^ + \A_\^-p^\A_\^} (B2) 



43 



where the modes A+, and A_ are related to the original perturbations by 
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1 - 
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^0. 



Pt 
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[b, sVgp't J J 

P . ^(V6 3gg (pi~3p|.)+8m^p|, ^\ ^ 



(B3) 



The modes A_|_, and A_ are the canonical variables of the system (they are the analogous of the Mukhanov- 
Sasaki [sl] variable v in the standard case of scalar field isotropic inflation). As in the standard case, their early 
times frequency is given by the momentum p, up to O (H/p) subdominant corrections. Since the momentum changes 
adiabatically at early times {p/p'^ — O (H/p)), we can set the initial conditions for the canonical modes according to 
the adiabatic vacuum prescription, precisely as done in the standard case [37| : 



H 



+ ■ 'in 



A_, 



1 



2p 



+ , in 



— ^-^,in — in — 



(B4) 



which are O {H/p) accurate. From eqs. (jB3p and (jB4p we thus obtain the initial conditions for {^P, a, di} and their 
time derivatives. Finally, the first, third, and fourth of eqs. (|56p provide the initial conditions for the nondynamical 
modes $, ^-nd do ■ In this way, we have the initial conditions for all the modes of the system ()59p . 
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